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ABSTRACT 

We present new HST/WFPC2 observations of 20 fields centered around T 
Tauri stars in the Chamaeleon I star forming region. Images have been obtained 
in the F631N ([OI]A6300), F656N (Ha) and F673N ([SII]AA6716, 6731) narrow- 
band filters, plus the Johnson V-band equivalent F547M filter. We detect 31 
T Tauri stars falling within our fields. We discuss the optical morphology of 
10 sources showing evidence of either binarity, circumstellar material, or mass 
loss. We supplement our photometry with a compilation of optical, infrared 
and sub-millimeter data from the literature, together with new sub-mm data 
for three objects, to build the Spectral Energy Distributions (SED) of 19 single 
sources. Using an SED model fitting tool, we self-consistently estimate a number 
of stellar and disk parameters, while mass accretion rates are directly derived 
from our Ha photometry. We find that bolometric luminosities derived from 
dereddened optical data tend to be underestimated in systems with high a 2 _24 
IR spectral index, suggesting that disks seen nearly edge-on may occasionally 
be interpreted as low luminosity (and therefore more evolved) sources. On the 
other hand, the same a 2 ~24 IR spectral index, a tracer of the amount of dust in 
the warmer layers of the circumstellar disks, and the mass accretion rate appear 
to decay with the isocronal stellar age, suggesting that the observed age spread 
(~ 0.5 — 5 Myr) within the cluster is real. Our sample contains a few outliers 
that may have dissipated their circumstellar disks on shorter time-scale. 

Subject headings: open clusters and associations: individual (Chamaeleon-I) - 
stars: pre-main-sequence — stars: luminosity function, mass function 



a Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 1-50125 Firenze, Italy 
department of Planetary Sciences, 1629 E. University Blvd., Tucson AZ 85721 
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Introduction 



The formation of stars, brown dwarfs, and planets is associated with a rich phenomenol- 
ogy spanning almost the entire electromagnetic spectrum. During the star formation process, 
the circumstellar disk material can be accreted into one (or more) central objects, ejected 
along the stellar polar axis, dispersed through photo-evaporation or condensed into plan- 
etesimals. Each of these processes is traced by a characteristic set of morphological and 



spectroscopic signatures, such as UV excess for stellar mass accretion (IHartmann et al. 



1998: Muzerolle et al. 


2001: 


Bouvier et a 




2007 


, shock-excited line- 


Herbig-Haro obiects ( 


Ballv et al. 


2000, 


2007; 


Reipurth et al. 


2003 


Podio et al. 2006). CO outflows from the disks 


Richer et al. 


2000: 



Whelan et al. 


2005; 


Arce et al. 


200 


7), or 



mid-IR and sub-mm signatures in the spectral energy di stribution due to grain growth and 



possible disk clearing by a forming planetary system fe.g. iRicci et al. 



2011 



Calvet et al. 



2002 



Dominik et al. 1120071 : lEspaillat et al 



2008 



2010 



Andrews et al. 



Blum fc Wurm M2008I : 



Pascucci fc Sterzik II2009I ). Each tracer provides critical information on different aspects of 
the protostellar system, but only a comprehensive view allows reconstructing the full evo- 
lutionary scenario. In particular, by combining optical and IR observations it is possible to 
address key open questions like the formation of stellar companions, the fraction of disk mass 



accreted into the central star vs. time, the times cale of planetesimal formation (IKrumholz 
20061 : iMever et al-lbood kpai & Lauretta I boioh . 



The Hubble and the Spitzer Space Telescopes offer the best available combination of 



field-o 
these i 


f-view, sensitivity, spatial resolu 


tion, and wave 


ength coverage for studying most of 


Dhenomena ( 


Padgett et al. 


1999; 


O'Dell & Wen 


1994J: 


Be 


Jlv et al. II200C 


; Robberto et al. 


2004; 


Ricci et al. 


200£ 


\\ Allen et al. < 


2004: Gutermuth et al. 


2004: Meeeat 


h et al. 2004: 


Muzerolle et al. 12004: . 


A.pai et al. 


200, c 


3; Pascucci et al. II2OO8 


>; 1 


^uhman I2OO8: 


Luhman et al. 



20041 : IColgan et al. 1 12007 



the a r chetype of star formi n g regions (lO'Dell &: Wen 111 994; B allv et al. Il2000t iRobberto et al. 



Ricci et al. 



2008; 



a and its associated young cluster, 



I lO'Dell fc Hennev IboOS ). HST observati 



cions 

resolved about 200 circumstellar disks, tens of jets and provided accurate broad-band pho- 
tometry needed to determine the fundamental stellar parameters. Unfortunately, the Orion 
cluster is relatively distant (~420pc), crowded and projected over the bright M42 HII region. 
Those factors make the Orion Nebula cluster a problematic target for Spitzer, while ground- 
based mid-IR observation s having adequate spa t ial resolution attain low sensitivity limits 



longward than ~ 3.5/j.m (IRobberto et al. I l2005t ISmith et al. II2005I ). On the other hand 



there are other closer regions that have been extensively investigated by Spitzer, but none of 
them has been studied with the HST with comparable detail. Among them, the Chamaeleon 
I region is possibly the best site for combining the unique Spitzer and HS T capabilities. It 
is one of the nearest star-forming regions (d=160-170 pc; for a review see iLuhman 1 120081 ) 
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and nearly coeval to Orion ( IDa Rio I l201dl . age ~2 Myr,), making its low- mass members 
approximately 10 times brighter than Orion. The cluster is young enough that it retains a 
significant population of primordial disks, but it is old enough that most of its members are 
no longer highly obscured by dust (typically A v < 2). Like Orion, because of the relatively 
low extincti on, optical wavelengths are accessible for the spectral classification of the stellar 
population ( Comeron et al. 12004 ; Luhman 2004, and references the rein) and for measuring 



accretion diagnostics (IMohanty et al. 



2005 



Muzerolle et al. 1 120051 ) . Optical and near-IR 



imaging and spectroscopic surveys of Chamaeleon I have produced an extensive and vir- 



tually complete census of the stellar and substellar cluster members in the field ( ILuhman 



20071 ): there are 237 known members, 33 of which have spectral types indicative of brown 
dwarfs (> M6). The IMF of Chamaeleon I reaches a maximum at a mass of 0.1-0.15 M , 
somewhat lower than the IMF peak in Orion. 

In this paper we report on a study aimed at probing a sample of Chamaeleon I sources 
with the Wide Field Planetary Camera 2 (WFPC2) onboard HST. Given the extent of 
region, about 5.5° x 1.5°, we targeted selected fields centered on brown dwarfs, class I and 
II Pre-Main-Sequence (PMS) objects measured by Spitzer to uncover substellar companions 
down to 15 AU separation and to directly image circumstellar disks and jets. We concentrate 
here on the sample of stellar PMS objects, leaving the discussion of the brown dwarf survey 
to a second paper (Luhman et al., in preparation). In Section 2 we illustrate our observing 
and data reduction strategy, detailing the extraction of photometry for both point sources 
and extended objects. In Section 3 we present the HST photometry, complemented by a 
compilation of IR and millimeter data available in the literature, including new data in 
the sub-mm range for 3 sources. We also list the main physical parameters of the stellar 
sources taken from the literature and derive mass accretion rates from our Ha photometry. 
In Section 4 we illustrate the morphology of individual objects, while in Section 5 we model 
the SEDs for 19 sources using a SED fitting tool; we derive disk parameters that can be 
compared wih stellar mass, age, luminosity and mass accretion rate. Finally, in Section 6 we 
summarize our findings. 



2. Observations 

2.1. HST Data Acquisition and Reduction 

The data presented in this paper have been obtained with the Wide Field Planetary 
Camera 2 (WFPC2) onboard the Hubble Space Telescope (HST) in early 2009 (HST GO 
program 11983, P.I. Robberto). These are among the latst data taken with the aging instru- 
ment, just preceeding the Servicing Mission 4 which replaced WFPC2 with WFC3. 
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We targeted 20 fields centered on T Tauri stars, detecting 18 of them0. For two sources, 
ISO- Chal 150 and ChaJl 108 1 938- 7731522, kn own for being highly obscured by circumstellar 



dust (jCambresy et al. I Il998l ; iLuhman 1 120071 ). we can only provide detection upper limits. 



Other 13 members of the Chamaeleon I complex lying in our imaged fields have also been 
identified and are presented in this paper. One HST orbit was dedicated to each field, placing 
the target on the standard aperture spot of the WFPC2 Planetary Camera (PC) chip, with 
45.5 mas/pixel scale. We did not constrain the telescope roll angle. 

Observations were carried out in the narrow-band filters F631N ([01] A6300) , F656N 
(Ha), and F673N (centered on the [SII] AA6716, 6731 doublet), plus the F547M medium- 
band filter roughly corresponding to the Johnson V-band. The single exposure times were set 
at 100 s (F631N), 40 s (F656N), 100 s (F673N) and typically 40 s for the F547M filter (only 
for the brightest stars the F547M were shorter to prevent saturation) with four exposures 
per filter taken in two groups of two. Each group was centered at a sli ghtly different position 



'two -point dither") for optimal bad pixels and cosmic ray rejection ( jBiretta fc McMaster 



2008|). 



Each set of four i mages was processed by the OPUS pipeline and combined using the 
MultiDrizzle software ( IFruchter &: Sosey Il2009l ). We used the MultiDrizzle parameters rec- 
ommended for two-point dithered observations, treating separately the PC chip from the 
other 3 Wide Field (WF) chips, due to the different pixel scale. 



2.2. Source Identification and Photometry 

The images processed by MultiDrizzle are corrected for geometric distortion introduced 
by the WFPC2 optics, cleaned from bad pixels and cosmic rays, and recombined into a single 
integrated image properly oriented in Right Ascension and Declination. We used STARFIND 
(in the STSDAS library of PYRAF) to determine the location on the CCDs of all sources in 
the field. 

After visually inspecting each individual source to reject false identifications we per- 
formed aperture photometry with DAOPHOT using a circular aperture of 3 pixels in radius 
(corresponding to 0'.'137 on the PC and 0'.'299 on the WF chips). The small extraction radius 
was chosen to optimally estimate the magnitude of weak sources, the large majority. We did 
not perform PSF photometry due to the non-negligible Charge Transfer Inefficiency trails of 



lr The HST archive contains four other fields, centered on 2MASSJ11095493-7635101, HN10E, IS0235 and 
IS079. They have been observed but the images appear compromised and have not been used in this work. 
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the aging instrument. Since the values of the zero points derived from the PHOTOFLAM 
keyword refer to counts measured within an "infinite" aperture^, we performed an aperture 
correction to convert the counts measured in our 3 pixel extraction radius to the value ex- 
pected in the 0'.'5 radius associated with the zero points. To this purpose we selected, for 
each filter and camera, a set of bright unsaturated and isolated sources and compared the 
results obtained with the two apertures. The average ratio, estimated with a sigma clip- 
ping algorithm, provided the aperture correction, together with the corresponding standard 
deviation. This allowed us to derive absolute magnitudes and errors in the HST STMAG 
system. For the narrowband filters we derived the flux in the more app ropriate PHOTFLAM 
(erg cm -2 s~ 1 A~ ;1 ") system (for the conversion from counts to flux see iBaggett et al. 1120021). 



Dolphin 



Finall y, we applied a correction for Charge Transfer Efficiency (CTE) loss following 
(ESS). 

A direct measure of the FWHM reveals that a number of sources are extended. To 
discriminate between extended and non-extended sources, we used the set of bona-fide point 
sources to build average Point Spead Functions valid for the PC and WF cameras. The 
FWHM of these PSFs turned out slightly larger than the theoretical FWHM provided by the 
HST PSF simulator Tiny Tim (of the order of 2 pixels for both channels), as expected since 
MultiDrizzle cannot fully recover the optimal PSF with only two pointings. The photometry 
of extended sources was determined using extraction radii large enough to contain all the 
signal above ~ 3a sky noise floor. We did not apply any corrections for CTE losses, as this 
is mostly relevant for small apertures. 



2.3. Sub-mm and mm observations 

As part of a larger project carried our with the Atacama Pathfinder Experiment (APEX) 
and with the Australia Telescope Compact Array (ATCA) facilities, we obtained new mil- 
limeter wavelength data for three sources: ESO Ha-559 (source #10) , ISO-Chal 10 (#5) and 
HH48 A (#7). The first two were observed with APEX, a 12 -m submm telescop e in Chile's 



Atacama desert, with the Large APEX BOlometer CAmera (jSiringo et al 1120091 LABOCA; 
) operating at the central frequency of 345 GHz (or 870 //m). The angular resolution of 
LABOCA is about 19 arcseconds and its total field of view is about 11 arcminutes. The ob- 
servations, conducted under the project 086.C-0653, were carried out on 2010 October 20-21 
for a total of 4.4 h (on source) for ESO Ha-559 and on 2010 October 29-30 for a total 4.8 



2 in the case of WFPC2 the flux in an infinite aperture is estimated to provide 1.096 times, i.e. a tenth of 
a magnitude, the counts measured in a 0'.'5 radius aperture. 
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h (on source) for ISO-Chal 10. Atmospheric conditions were excellent, with a precipitable 
water vapor around 0.34 mm, corresponding to a zenith opacity of 0.18. The sky opacity was 
measured every hour with skydips. The pointing of the telescope was checked every hour on 
the nearby quasar PKS1057-79. The absolute flux calibration was performed by observing 
the secondary calibrator B13134 every one or two hours. The telescope focus was checked by 
observing the star r\ Carinae, at least once per day. The observations were performed on the 
fly mode with a rectangular scanning pattern. The data were re duced with the BOl ometer 



Array Analysis Software, following the procedures described in ISiringo et al I (120091 ). The 
total (flux) calibration error is about 20% and the rms noise level in the final maps was 
~5 mJy We report a clear detection of the submm emission from ESO Ha-559, but only an 
upper limit could be derived from our data for ISO-Chal 10. 

HH 48 A was observed at 3.3 mm with ATCA and the new CABB digital filter bank, 
which provides a total continuum bandwidth of 4 GHz. Observations were carried out at a 
central frequency of 91.000 GHz (3.294 mm) on 2009 Oct 14. The ATCA array was in the 
hybrid H168 configuration, providing an angular resolution of about 3 arcsec at 3.3 mm. The 
gain was calibrated with frequent observations of 1057-797. The passband was calibrated 
using 1921-293, and the absolute flux scale was determined through observations of Uranus. 
The uncertainty on the ATCA calibrated flux is about 30% at 3.3 mm. The MIRIAD package 
was used for visibilities calibration, Fourier inversion, deconvolution and imaging. The rms 
noise on the ATCA map was about 0.45 mJy. 

We also searched the literature finding seven more sources previously detected at wave- 
lengths longer than 100 fim. All fluxes are reported in Table HJ illustrated in the next section. 



3. Results 

3.1. Final source list and HST photometry 

The location of our 20 WFPC2 fields, each composed of 3 WF images of ~ 72"8 x 72 / .'8 
and 1 PC image of ~ 36'.'4 x 36'.'4, is shown in Figure [TJ A few fields in the northern region 
appear partially superimposed, but they have been processed and analyzed separately. To 
facilitate data retrieval, in Tabled] we provide, for each field, the original HST visit number 
(by order of execution) and target name reported by the HST data archive, together with 
the more common SIMBAD target name used in this paper. 



Besides our 18 detected ta rgets, other 13 object s previously known to be members of 
the Chamaeleon I association (jLuhman et al. 1 120081 ) ended up in our imaged fields. For 
each of these 31 sources, Table [2] provides: an entry number (column 1), the coordinate- 
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Table 1. Main target list 



Visit 


HST Target 


Source name 


17 


T14A 


HH 48 IRS 


18 


IS0225 


ISO-Chal 225 


19 


OTS32 


ISO-Chal 232 


20 


T14 


CT Cha 


21 


IRN 


ISO-Chal 150 a 


22 


T47 


HBC 584 


23 


CHSM15991 


CHSM15991 


24 


CRHF574 


ESO-Ha 574 


25 


CRHF569 


ESO-Ha 569 


26 


CHAJ11081938- 7731522 


Cha J11081938-7731522 a 


27 


2MASSJ10533978- 7712338 


2MASS J10533978- 7712338 


28 


HN21E2 


Hn 21E 


29 


T12 


ISO-Chal 10 


30 


T42 


CED 112 IRS 4 


31 


T3A 


SX Cha 


32 


T5 


Ass Cha T 2-5 


33 


IS0252 


ISO-Chal 252 


34 


CHXR20 


UX Cha 


35 


T16 


Ass Cha T 2-16 


36 


CRHF559 


ESO-Ha 559 



a Main target not detected. 
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Fig. 1. — All WFPC2 fields observed in our progr am superimposed to a map of the A 



extinction in Chamaeleon I (jCambresy et al. 1119971 ). The highest contour corresponds to 
Aj = 10.7 while the lowest contour to Aj = 0. The label of each field provides the visit 
number and a suffix, either TT (T Tauri star) or BD (Brown Dwarf), to indicate the nature 
of the primary target. In this paper we concentrate on the TT fields. The random orientation 
of the characteristic WFPC2 chevron pattern is due to unconstrained HST roll angle. For 
the association between HST visit and the corresponding target see Table 1. 
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based 2MASS source identifier, when available (column 2); an alternate ID taken from the 
literature (column 3); equatorial coordinates from the original fits header (column 4 and 5); 
the HST visit, and the (x,y) coordinates on the drizzled fits files (column 6 and 7). The 
accuracy of the celestial coordinates, driven by the absolute positions of the HST guide stars, 
is typically about 0'/25. 



Table 2. Observed members of the Chamaeleon I association 



Nr. 2MASS SIMBAD R.A. (J2000.0) Decl. (J2000.0) Visit FITS COORD. (X-Y) 



1 


2MASS 


J10533978- 7712338 




10 53 39.78 


-77 12 33.9 


27 


1674-2119 


2 






Ass Cha T 2-3 B 


10 55 59.73 


-77 24 39.9 


31 


1883-2291 


3 






SX Cha 


10 55 59.76 


-77 24 40.1 


31 


1837-2245 


4 


2MASS 


J10574219-7659356 


Ass Cha T 2-5 


10 57 42.20 


-76 59 35.7 


32 


1705-2153 


5 


2MASS 


Jl 1025504-7721508 


ISO-Chal 10 


11 02 55.05 


-77 21 50.8 


29 


1681-2143 


6 


2MASS 


Jl 1040909-7627193 


CT Cha 


11 04 09.09 


-76 27 19.4 


20 


1631-2084 


7 


2MASS 


J11042275- 7718080 


HH 48 A 


11 04 22.69 


-77 18 09.1 


17 


1545-1991 


8 






HH 48 B 


11 04 23.31 


-77 18 07.5 


17 


1502-2024 


9 


2MASS 


Jl 1045701-7715569 


Ass Cha T 2-16 


11 04 57.01 


-77 15 56.9 


35 


2346-2647| 


10 


2MASS 


Jl 1062554-7633418 


ESO-Ha 559 


11 06 25.55 


-76 33 41.9 


36 


2382-265£ 


11 


2MASS 


Jl 1064510-7727023 


UX Cha 


11 06 45.10 


-77 27 02.3 


34 


2381-2665 


12 


2MASS 


Jl 1081648-7744371 


Ass Cha T 2-34 


11 08 16.49 


-77 44 37.2 


21 


1551-577 


13 


2MASS 


J11081703- 7744118 


ISO-Chal 137 


11 08 17.03 


-77 44 11.8 


21 


1534-832 


14 






Cha J11081938- 7731522 


11 08 19.38 


-77 31 52.2 


26 


not detected 


15 






ISO-Chal 150 


11 08 37.1 


-77 43 51 


21 


not detected 


16 


2MASS 


Jl 1094525-7740332 


ISO-Chal 201 


11 09 45.26 


-77 40 33.3 


23 


983-972 


17 


2MASS 


Jl 1094621-7634463 


Hn 10E 


11 09 46.21 


-76 34 46.4 


30 


1349-996 


18 


2MASS 


Jl 1095262-7740348 


CHSM 15991 


11 09 52.62 


-77 40 34.9 


23 


1634-2091 


19 


2MASS 


Jl 1095340-7634255 


CED 112 IRS 4 


11 09 53.41 


-76 34 25.5 


30 


2401-2638 


20 


2MASS 


J11095437-7631113 


ISO-Chal 225 


11 09 54.38 


-76 31 11.4 


18 


1585-2054 


21 


2MASS 


J11100010-7634578 


WW Cha 


11 10 00.11 


-76 34 57.9 


19 


1212-329 


22 


2MASS 


J11100369-7633291 


ISO-Chal 232 


11 10 03.69 


-76 33 29.2 


19 


2386-2684 


23 


2MASS 


J11100469-7635452 


Cha T 2-45a 


11 10 04.69 


-76 35 45.3 


30 


704-405 
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In Table [3] we provide the WFPC2 photometry together with the photometric extraction 
area for extended sources. If a source is not detected we report the 3cx upper limit estimated 
over a 3 pixel aperture radius. 



Table 2 — Continued 



Nr. 




2MASS 




SIMBAD 


R.A. (J2000.0) 


Decl. (J2000.0) 


Visit 


FITS COORD. (X-Y) 


24 


2MASS 


Jll 102852- 


7716596 


Hn 12W 


11 10 28.52 


-77 16 59.6 


22 


1366-1401 


25 


2MASS 


J11103481- 


-7722053 


[LES2004] a Chal 405 


11 10 34.81 


-77 22 05.3 


33 


1274-426 


26 


2MASS 


J11104141- 


-7720480 


ISO-Chal 252 


11 10 41.42 


-77 20 48.1 


33 


2313-2633 


27 


2MASS 


Jll 104959- 


7717517 


HBC 584 


11 10 50.00 


-77 17 51.8 


22 


1463-1921 


28 


2MASS 


J11105076- 


7718031 


ESO-Ha 568 


11 10 50.77 


-77 18 03.2 


22 


630-763 


29 


2MASS 


J11111083- 


-7641574 


ESO-Ha 569 


11 11 10.83 


-76 41 57.4 


25 


1494-1950 


30 


2MASS 


Jl 1142454- 


-7733062 


Hn 21W 


11 14 24.54 


-77 33 06.2 


28 


2441-2862 


31 


2MASS 


J11142611- 


-7733042 


Hn 21E 


11 14 26.11 


-77 33 04.3 


28 


2328-2905 


32 


2MASS 


J11145031- 


-7733390 


BYB 53 


11 14 50.32 


-77 33 39.0 


28 


278-980 | 


33 


2MASS 


J11160287- 


-7624533 


ESO-Ha 574 


11 16 02.88 


-76 24 53.3 


24 


1558-2022 g 



Lopez-Marti et al. I (12004f l 



Table 3. Photometry with F547M, F631N, F656N and F673N filters 



IN . 


r o^t i iv± 
[mag] 


r DO ±1N 

[iu eig s em /\ j 


[iu eig s em i\ 


r u i oin 

If! -16 orcr nrn~ 2 

iu eig s em /\ 


iNoies 


1 


19.38 ± 0.18 


-0.37 


1.6 ± 0.9 


1.1 ± 0.5 




2 


16.14 ± 0.08 


11 ± 3 


37 ± 9 


16 ± 4 




3 


16.71 ± 0.03 


15 ± 0.8 


66 ± 4 


15.6 ± 0.8 


extended; 0'.'5 aperture 


4 


14.80 ± 0.03 


31.2 ± 1.7 


195 ± 10 


40 ± 2 


binary: / /5 aperture 


5 


16.25 ± 0.08 


10 ± 3 


62 ± 14 


12 ± 3 




6 


13.00 ± 0.07 


>230 


>770 


>250 


companion; 0'.'5 aperture 


7 


17.89 ± 0.03 


4.1 ± 0.3 


17.4 ± 1.1 


3.8 ± 0.2 


disk; 0'.'5 aperture 


8 


21.3 ± 0.3 


0.32 ± 0.16 


2.8 ± 0.8 


0.44 ± 0.11 


disk; 1" aperture i 


9 


18.91 ± 0.04 


1.75 ± 0.15 


3.4 ± 0.4 


2.27 ± 0.16 


extended; 0'/5 apertiire 


10 


20.7 ± 0.3 


-0.35 


-1.3 


-0.43 


l 


11 


14.78 ± 0.07 


88 ± 18 


130 ± 30 


100 ± 20 




12 


16.53 ± 0.05 


11 ± 3 


27 ± 5 


15 ± 3 




13 


18.51 ± 0.11 


1.7 ± 0.7 


5.3 ± 1.7 


2.4 ± 0.8 




14 


>24 


<0 035 


<0 14 


<0 023 




15 


>24 


<0.037 


<0.15 


<0.025 




16 


19.94 ± 0.02 


0.62 ± 0.06 


2.4 ± 0.2 


0.86 ± 0.06 




17 


18.41 ± 0.10 


2.6 ± 0.9 


20 ± 4 


3.4 ± 1.0 




18 


21.47 ± 0.09 


0.20 ± 0.04 


0.32 ± 0.13 


0.24 ± 0.03 




19 


19.50 ± 0.11 


2.9 ± 0.2 


9.85 ± 0.7 


2.15 ± 0.15 


extended; 0'/5 aperture 


20 


21.6 ± 0.5 


-0.13 


-0.33 


-0.17 




21 


14.10 ± 0.02 


148 ± 8 


650 ± 30 


1.49 ± 8 


extended; 1" aperture 


22 


18.78 ± 0.14 


2.6 ± 1.0 


4.6 ± 1.8 


4.2 ± 1.3 
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Table |U provides a compilation of near-IR photometric data for the 31 detected sou rces, 
either from the ground or from Spitzer ( ILuhman et al. 1120081 ; iLuhman fc Muench 1120081) . We 
also li st fluxes at 0.87 mm, 1.3 mm and 3.3 mm from previously published data (Bello che et al. 



(2011), using APEX/LABOCA at 870 micron and ~ 19" angular resolution; iHenning et al. 



19931 ). using SEST at 1.3 mm and ~ 23" angular resolution), together with our new data 



for 3 sources (see Section 3.2). The appropriate reference for each individual source is given 
in the Appendix. 

In Table Owe present a compilation of the main physical parameters of the stellar sources 
(spectral type, effective temperature, extinction and bolometric luminosity) reported in the 
literature, together with the equivalent width of the Ha line, accretion luminosity and mass 
accretion rate, estimated as follows. We have computed for each star the photospheric 
continuum in the F656N filter and subtracted it from the measured flux. The continuum 
ha s been evaluat e d per forming synthetic photometry. The set of BT- Settle synthetic spectra 
of lAllard et al. I (12010|) was fi r st inte rpolated at the T e g of each source and then reddened 
using the Aj from ILuhman I (120071 ) . For each spectrum we computed the photospheric 
(F547M — F656iV) color; then using the measured F547M magnitude as a reference, we 
rescaled the photometry to observed fluxes to derive the photospheric continuum F656N . 
The Equivalent Width (E.W.) of the H a excess was then derived from the ratio between 
the flux excess and the continuum, multiplied by the equivalent width of the F656N filter 
profile. We converted the H a excess in units of stellar luminosity by estimating the fraction 
of stellar (photospheric) flux entering in the F656N filter window, using once again the 
BT-Settle spec tra. Then, knowing the bolometric luminosity of the sources from (from 
Luhman 1120071 see Table 4), we derived the H a excess in units of solar luminosity L & . The 
H a luminosit y can be then related to the overall accretion luminosity L accr assuming the 



formula from |Pe Marchi et al. I ( 120101 ) : 



logL accr = 1.72 + \ogL Hc 



and therefore the mass accretion rates using the relationship IGullbring et al. I ( 119981 ) 



0.8 ■ 



GM*M acc 



(1) 



(2) 



under the assumption that the mass infall onto the stellar surface starts from a distance of 
about 5 stellar radii. 



Table 3 — Continued 



N. 


F547M 
[magj 


F631N 

r-i n — 16 — 1 —2 X —11 

[10 erg s cm * A 1 J 


F656N 

r-i n — 16 —1 —2 X —11 

[10 D erg s cm z A i J 


F673N 

1 n — 16 —1 —2 X —11 

10 D erg s cm z A X J 


Notes 


23 


14 68 ± 03 


70 ± 12 


94 ± 17 


81 ± 13 




24 


1 7 Q3 + OQ 


2 8 + 10 


8 + 2 


3^ + 10 




25 


20 9 ± 3 


—0 34 


—1.2 


—0.52 




26 


-22.4 


-0.094 


2.4 ± 1.2 


-0.13 




27 


16.64 ± 0.08 


10 ± 3 


49 ± 11 


13 + 3 




28 


18.27 ± 0.03 


2.72 ± 0.18 


6.3 ± 0.5 


3.4 ± 0.2 




29 


20.28 ± 0.14 


0.21 ± 0.15 


2.4 ± 0.7 


0.35 ± 0.12 


disk; 1" aperture 


30 


18.20 ± 0.12 


2.5 ± 1.0 


25 + 6 


3.8 ± 1.2 




31 


20.7 ± 0.3 


-0.28 


-0.88 


-0.34 




32 


<15 


67 ± 4 


99 + 5 


78 + 4 




33 


19.60 ± 0.07 


1.28 ± 0.18 


2.7 ± 0.6 


1.96 ± 0.18 


disk; 1" aperture 



Table 4. IR and sub-mm photometry 



N. 


R a 
SED 












J° 


H c 


Ks c 


3.6 i-im 


4.5 fim 


5.8 fim 


8.0 fim 


24 


fATXl 


F0.87mm 




F 3 .3mrn 




[mag] 




[mag 


i 




[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mJy] 


[mJy] 


[mJy] 


1 




15 


.56 


± 





06 


1 O OO _i_ A AO 

lo.zo ± 0.03 


1 O 1 A _l_ A AO 

LZ. 14 zh U.IM 


11 CO _J_ A AO 

11. do zh V.VZ 


11 ACl _l_ A AO 

11.49 zh V.UZ 


1 1 An _L a ao 
11. uy ± 0.02 


1 n OA _l_ A AO 

1U.DU ± U.U3 


A TO _l_ A A A 
9. (O ± 0.04 


5.28 


± 0.04 








2 




12 


.35 


± 





02 


























3 




12 


.35 


± 





02 


1 a a c: _i_ a ao 


A O A _L A AO 


O Cfl _l_ A AO 

o.o9 zh \).\)Z 




















4 




12 


17 


± 





0.3 


1 A AO _1_ A AO 

lU.4d zh U.Uz 


A EC _L_ A AO 


A OC _J-_ A AO 

y.zo zh iiAiZ 


O *7A _l_ A AO 

o. / 4 zh V.VZ 




O A V A AO 

0.0 ( ±0.03 




4.35 


± 0.04 








5 


15. 1 < ± 0.05 


13 


.40 


± 





05 


11 C ^ _l_ A AO 

11. do ± 0.03 


1 A O/? _l_ A AO 

10. oo zh U.Ud 


1 A Ad _i_ A AO 

10.45 zh \i.\iZ 


A TO _l_ A AO 

9. io zh V.VZ 




A OA -L- A AO 

y.20 ± o.oo 


o a C _1_ A AO 

O.OO ± O.Oo 


5.90 


± 0.04 


< 14.8 






(i 














A T1 _l_ A AO 

y.7i ± u.uz 


O A _J_ A A K 

o.y4 ± 0.05 


O _]_ A AO 

o.OD zh U.UZ 










2.44 


± 0.04 








7 

8 














16. 42^ 
18. 43 4 


14. 32 4 
15. 66 4 


12. 54 4 

13. 85 4 


10.06 zt 0.02 


9.23 1 0.02 


8.25 ± 0.03 


6.97 ± 0.04 


3.47 


± 0.04 






3.95 


9 


16.65 0.05 


14 


.68 


± 


0. 


04 


12.17 ± 0.02 


10.97 zt 0.02 


10.41 zt 0.02 


9.80 zh 0.02 


9.57 ± 0.02 


9.18 ± 0.04 


8.63 ± 0.04 


6.11 


± 0.04 








10 




15 


.88 


± 


0. 


00 


13.01 zt 0.03 


12.01 zh 0.02 


11.49 zt 0.02 


10.83 zt 0.02 


10.42 ± 0.02 


10.07 ± 0.03 


9.38 ± 0.04 


5.34 


± 0.04 


44.1 






11 




12 


07 


± 


0. 


02 


10.18 zt 0.02 


9.20 zt 0.02 


8.880 zt 0.019 


8.51 zt 0.02 


8.36 ± 0.02 


7.93 ± 0.03 


6.94 ± 0.04 


4.45 


± 0.04 








12 




13 


12 


± 


0. 


03 


11.20 ± 0.03 


10.34 zh 0.03 


10.02 zt 0.02 


9.78 zt 0.02 


9.75 ± 0.02 


9.73 ± 0.03 


9.71 ± 0.03 












13 




14 


34 


± 


0. 


04 


11.79 ± 0.03 


11.06 zt 0.03 


10.67 zt 0.02 


10.25 zt 0.02 


10.18 ± 0.02 


10.13 ± 0.03 


10.15 ± 0.04 












16 


18.02 ± 0.05 


15 


.37 


± 


0. 


00 


12.36 zt 0.03 


11.45 zh 0.02 


11.03 zt 0.02 


10.55 zt 0.02 


10.40 ± 0.02 


10.39 ± 0.03 


10.39 ± 0.04 


9.82 


± 0.18 








17 




14 


72 


± 


0. 


05 


11.95 zt 0.02 


10.74 zh 0.02 


10.05 zt 0.02 


9.49 zt 0.02 


8.92 ± 0.02 


8.46 ± 0.03 


7.59 ± 0.04 


3.70 


± 0.06 


1552 






18 














16.05 ± 0.11 


14.87 zh 0.07 


14.13 zt 0.07 


11.98 zt 0.02 


11.37 ± 0.02 


10.71 ± 0.03 


9.87 ± 0.03 


7.21 


± 0.05 








19 


15.10 ± 0.05 


13 


.52 


± 


0. 


05 


9.47 =t 0.02 


7.79 ± 0.05 


6.46 zt 0.03 






3.72 ± 0.03 


3.04 ± 0.03 






71 






20 




17 


27 


± 


0. 


12 


15.05 =t 0.11 


13.80 zh 0.10 


13.14 =t 0.07 


11.22 zt 0.02 


10.35 ± 0.02 


9.61 ± 0.07 


8.70 ± 0.04 


4.85 


± 0.04 


<290 






21 




10 


.95 


± 


0. 


04 


8.71 =t 0.03 


7.21 ± 0.08 


6.08 zt 0.05 






3.70 ± 0.03 


3.15 ± 0.04 






1501 


407.9 l— J 




22 




14 


72 


± 


0. 


05 


11.77 =t 0.03 


10.26 zt 0.02 


9.44 zt 0.02 


8.18 zt 0.02 


7.52 ± 0.02 


6.84 ± 0.03 


6.04 ± 0.04 


3.23 


± 0.04 


<290 




23 




12 


.39 


± 


0. 


04 


10.56 ± 0.02 


9.64 =t 0.02 


9.24 zt 0.02 


8.70 zt 0.02 


8.41 ± 0.02 


8.10 ± 0.03 


7.55 ± 0.03 


4.58 


± 0.04 








24 




13 


.97 


± 


0. 


05 


11.73 ± 0.02 


11.11 zt 0.02 


10.78 zt 0.03 


10.37 =t 0.02 


10.23 ± 0.02 


10.26 ± 0.03 


10.26 ± 0.04 


10.00 


± 0.09 








25 




15 


(il 


± 


0. 


06 


12.04 ± 0.02 


10.72 zt 0.02 


10.03 zt 0.02 


9.60 zt 0.02 


9.43 ± 0.02 


9.42 ± 0.03 


9.42 ± 0.04 


9.19 


± 0.10 








26 


20.06 ± 0.05 


17 


27 


± 


0. 


05 


13.86 ± 0.03 


12.89 zt 0.03 


12.27 zt 0.02 


11.45 zt 0.02 


11.02 ± 0.02 


10.57 ± 0.03 


9.75 ± 0.03 


7.05 


± 0.04 








27 




13 


84 


± 


0. 


05 


11.15 ± 0.02 


9.95 zt 0.02 


9.17 zt 0.02 




8.13 ± 0.02 




6.37 ± 0.04 


2.96 


± 0.04 


44.95 






28 


16.65 ± 0.05 


14 


.50 


± 


0. 


05 


12.04 ± 0.02 


11.10 zt 0.02 


10.75 zt 0.02 




10.26 ± 0.02 




10.22 ± 0.05 












29 




17 


34 


± 





12 


15.95 ± 0.09 


15.06 zt 0.09 


14.58 zt 0.10 


14.21 =t 0.03 


13.76 ± 0.03 


13.23 ± 0.06 


12.50 ± 0.05 


7.12 


± 0.04 


72 






30 




14 


24 


± 


0. 


04 


11.98 ± 0.05 


11.09 zt 0.05 


10.65 zt 0.04 


10.06 zt 0.04 


9.84 ± 0.04 


9.54 ± 0.04 


8.96 ± 0.04 


6.42 


± 0.04 








31 




15 


.47 


± 


0. 


05 


12.76 zt 0.03 


11.97 =t 0.03 


11.49 zt 0.02 


11.01 zt 0.04 


10.89 ± 0.04 


10.79 ± 0.06 


10.57 ± 0.08 












32 




11 


.83 


± 


0. 


05 


10.48 zt 0.03 


9.75 zt 0.02 


9.55 zt 0.02 


9.35 zt 0.02 


9.28 ± 0.02 


9.26 ± 0.03 


9.24 ± 0.03 


9.06 


± 0.05 








33 




17 


.25 


± 


0. 


12 


15.80 ± 0.07 


14.97 =t 0.08 


14.61 zt 0.11 





















lLopez-Marti et al~l i 2004ft 

' lLopez-Marti et al~l I2004D ; Second DENIS Release. 

C 2MASS Point Source Catalog. 

lHaiach et al. I i2004h , labeled as Cha I T14a 



Note. — Data at A > 3.6fim arc from Luhman et al. 12008): Luhman & Mucnch 12008); Bclloche et al. 12011): Hcnning ct al. 11993) ; this paper. 
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Table 5. Source Physical Parameters 



N. 


Sp. Type a 


T e // a 


Aj a 


Lf>oZ a 


Ha excess E.W. 


log L acc 


log M acc 






[K] 


[mag] 


M 


[A] 


[L S ] 


[M Q /yr] 


1 


M2.75 


3451 


0.63 


0.032 


11.9±23.5 


-3.14±0.47 


-10.3±0.5 


2 


M3.5 


3342 


0.79 


0.42 


8.5±9.4 


-2.23±0.32 


-8.8±0.3 


3 


MO 


3850 


0.79 


0.42 


122.8±10.8 


-0.85±0.04 


-7.8±0.1 


4 


M3.25 


3379 


0.34 


0.33 


47.5±5.8 


-1.58±0.05 


-8.2±0.1 


5 


M4.5 


3198 





0.081 


75.0±24.6 


-2.04±0.12 


-8.8±0.1 


6 


K5 


4350 


0.45 


0.95 








7 


K7 


4060 


0.45 


0.013 


124.4±10.8 


-2.23±0.04 


-10.0±0.1 


9 


M3 


3415 


1.6 


0.21 


0.9±4.3 


-3.49±0.77 


-10.3±0.8 


10 


M5.25 


3091 


1.01 


0.052 


39.6±26.1 


-2.56±0.22 


-9.3±0.2 


11 


K6 


4205 


1.08 


1.1 


19.2±11.2 


-1.05±0.20 


-8.1±0.2 


12 


M3.75 


3306 


0.29 


0.15 


22.3±9.7 


-2.28±0.16 


-9.0±0.2 


13 


M5.5 


3058 


0.32 


0.083 


24.3±17.4 


-2.63±0.23 


-9.2±0.2 


16 


M5.75 


3024 


0.5 


0.058 


52.0±7.4 


-2.44±0.06 


-9.0±0.1 


17 


M3.25 


3379 


1.01 


0.15 


140.6±38.9 


-1.39±0.11 


-8.3±0.1 


18 


M3 


3415 


0.79 


0.0029 


21.7±20.8 


-3.92±0.29 


-11.6±0.3 


19 


K5 


4350 


1.47 


3 


202.1±27.6 


0.44±0.05 


-6.5±0.1 


90 
zu 


Ml 7^ 


ooyo 


1 OA 


n m ^ 

U.UIO 


iy. i ztzo.u 


q oQ-i_n ^fi 

-o.zyzcu.oo 


i n s-i-n a 


21 


K5 


4350 


1.35 


5.5 


91.4±6.0 


0.42±0.03 


-6.6±0.1 


22 


K8 


3955 


2.14 


0.66 


13.7±17.0 


-1.30±0.35 


-8.7±0.3 


23 


Ml 


3705 


0.54 


0.43 


10.2±7.0 


-1.91±0.23 


-8.9±0.2 


24 


M5.5 


3058 





0.066 


28.5±14.8 


-2.65±0.18 


-9.2±0.2 


25 


M4 


3270 


1.91 


0.32 


19.5±19.4 


-2.02±0.30 


-8.6±0.3 


2G 


M6 


2990 


0.97 


0.022 


599.2±369.4 


-1.73±0.21 


-8.6±0.2 


27 


M2 


3560 


1.17 


0.42 


49.0±18.3 


-1.38±0.14 


-8.2±0.1 


28 


M4.25 


3234 


0.81 


0.11 


12.9±3.8 


-2.68±0.11 


-9.4±0.1 


29 


M2.5 


3488 


0.68 


0.003 


98.8±40.9 


-3.28±0.15 


-10.9±0.2 


30 


M4 


3270 


0.72 


0.11 


136.5±43.6 


-1.64±0.12 


-8.4±0.1 


31 


M5.75 


3024 


0.56 


0.042 


34.8±26.0 


-2.66±0.24 


-9.5±0.2 


32 


M2.75 


3451 





0.24 








33 


K8 


3955 


0.45 


0.0034 


82.9±25.5 


-3.04±0.12 


-10.99±0.11 



iLuhman I (|2007l ): iLuhman fc Muenchl |2008l) 
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4. Source Morphology as Revealed by HST Data 

In this section we illustrate our findings on the only 10 sources that appear either 
extended, binaries, or associated with jets and Herbig-Haro objects in the immediate vicinity. 
The other 23 sources appear point-like in our images. Our dataset is not homogeneous 
because of two main factors: a) the different sampling scales of the WFPC2 pixels (about 
50 mas/pixel for the PC chip, where the main targets were located, vs. 100 mas/pixel for 
the WF chips where most of the other sources are found); b) the presence of CTE tail, 
which depends on the position and brightness of the source on the chip. To mitigate the risk 
of misinterpreting extended features, we have indicated the direction of the CTE deferred- 
charge trails with an arrow in each image. For the brightest sources, i.e. Ass Cha T2-16 
(#9), CED 112 IRS4 (#19), and WW Cha (#21), we have also plotted on each image 
the isophotal contours of the corresponding PSF, for an immediate comparison. All images 
shown have the standard orientation with North up and East to the left. 



4.0.1. SX Cha (#3) 



The T Tauri star SX Cha, spectral type MO. 5 (ILawson. Feigelson and Huenemoerder 



19961 ). has a companion of comparable brightness ~2'.'l at P.A.~310° (INatta et al. 



200 



3). 



Spitzer data show the characteristic amorphous silicate emission of circumstellar disks at 10 
and 20 /xm, whereas the spectral features ob served in the 33-35 /xm rang e are characteristic 



of crystalline enstatite and forsterite grains (IKessler-Silacci et al. II2006I ). Our HST images 



(Fig. [2]) show SX Cha (source on the left) with its companion at ~2'/2 (350 AU) distance. 
The protuberance on the western side of SX Cha visible in the F673N ([SII]) image with a 
length of ~0V55 (87 AU), cannot be attributed to the CTE losses. The line intensity, together 
with the narrow and twisted morphology, suggests that it is a collimated jet. In this case, 
the circumstellar disk surrounding SX Cha would be oriented somewhat perpendicularly to 
the plane of the binary orbit. 



4.0.2. Ass Cha T 2-5 (#4) 



Detecte d for the first tim e by ISchwartz I (119771 ). Ass Cha T 2-5 has spectral type M3.25 
according to iLuhman I (120071 ). Our images (Fig. [3]) resolve this source in a close binary with 
separation ~0'/15, corresponding to ~25 AU projected distance at 160pc. The southern star 
appears brighter in all filters with the exception of Ha, where the northern source strongly 
dominates. This may indicate that the northern component has been observed in a phase 
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Fig. 2. — SX Cha observed in the PC chip (45.5 mas/pixel scale). Image centered at 
R.A.=10:55:59.53, Decl.=-77:24:40.3 (J2000.0). 



of strong accretion activity. The ratios between the peak counts of the northern vs. the 
southern star are approximately 0.57 (V-band), 0.89 ([01]), 2.8 (Ha) and 1.1 ([SII]). 



4.0.3. CT Cha (#6) 



CT C ha has been initially class ified as an emission-line star because of its Ha line 
variations (IHenize &: Mendoza 1119731) and later as a classical T Tauri star on the basis of 
its strong IRAS excess (IGauvin fe Strom II19921). ISO data have shown evidence of silicate 
emission in a circumstellar disk ( INatta et al. 1120001). The va riations of Ha line emission have 



been interpreted as accretion si gnatures by lHartmann et al. I ( 119981 ) . compatible with early 
observations of spectral veiling ( Rydgren 198ol ). More recent o bservations of CT Ch a have 
revealed the presence of two faint companions close to the star ( jSchmidt et al. 1120081 ). The 
first (source CI), about 6.3 mag fainter than the primary in the Ks-band, is located 2'/670 
± 0'.'038 at P.A.~315° of the star, corresponding to ~430 AU projected dis t ance at 160 pc. 
The second (C2) is about 2" at P.A.~45° of the primary. ISchmidt et al. I (120081 ) classified 
CI as physically associated with the primary due to its common proper motion, whereas 
for C2 they concluded that it must be a background object. In our PC images (Fig. H|) 
source CI is not detected while C2 is clearly visible at l'. / 96±0'. / 05 projected distance (Fig. 
13 right). We have found on the ESO data archive an H-band image taken with NACO in 
February 2006 in which only CI is visible (Fig. El left). A second image taken two years 
later, also with NACO but in the Ks filter, shows both CI and C2, with C2 at a distance 
of l / /91±0 / /03 from the primary (Fig. [5j center). Source C2 in our images is visible only in 
the [OI] and [SII] filters (Fig. [5j right). This, together with the fact that C2 (like source 
CI) has not changed position seems to indicate that it is indeed physically associated with 
V* CT Cha. Its nature, however, remains enigmatic. Finally, the point-like source in the 
[OI] image about l'/5 to the south of the main source may be real, as it cannot be easily 
attributed to a filter ghost. 



4.O.4. EE 48 A and B (#7, #t 



The classic Herbig-Haro object HH 48 (ISchwartz 1119770 is composed b y two close con- 
densations, designated as HH 48 A and HH 48 B. IWang &: Henning I (120061 ) found that HH 
48 A and B are elongated in directions that are roughly perpendicular, while a set of newly 
detected features (labeled C to F) are aligned along a direction pointing to HH 48 A. This 
sugg ests the presence o f two outflows driven from an embedded source, probably a binary 
star jBallv et al. Il2006[ ). Our images (Fig. [6J show only HH 48 A and B, i.e. HH 48 C-F are 
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Fig. 3. — Ass Cha T 2-5 observed in the PC chip. The countours levels are given 0.75, 
0.5 and 0.25 of the peak counts. Image centered at R.A. = 10:57:42.22, Decl.=- 76:59:36.4 
(J2000.0). 
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Fig. 4.— CT Cha observed in the PC chip. Image centered at R.A.=11:04:09.00, Decl.= 
76:27:19.5 (J2000.0). 
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not detected. In all niters, HH48 A appears as a bright point source associated with faint 
extended emission. While the elongation on the East side may be contaminated by CTE 
losses, the elongation on the West side is unambiguous. HH 48 B, located at ~2'.'5 (400 
AU) to the northeast of HH 48 A, is fainter and elongated nearly east- west, with a length of 
~0.7", in the V-band filter. In the [01] and Ha filters, however, it appears nearly unresolved. 
We speculate that the morphology of HH 48 B is compatible with the presence of a disk seen 
nearly edge-on, as the V-band elongation could be attributed to scattered light from a disk 
face, whereas the fainter line emission could arise from the inner disk region. Overall, our 
HST images confirm that HH 48 A and B point to different directions, forming an angle 
close to -30°. 



4.0.5. Ass Cha T 2-16 (#9) 



This source h as be en classified as a M3 emission-line star with a mass of ~0.26 M Q 
20081 ). Our images (Fig. [7]) indicate that Ass Cha T 2-16 has been 



flLafreniere et al. 



resolved by the HST, in particular in the [OI] where the FWHM~4 pixels is twice the size 
of the FWHM of unresolved sources. The SED of this source is shown in Section [5j 



4.0.6. CED 112 IRS 4 (#19) 



CED 112 IRS 4 is a T Tauri sta r located in a region rich of HH objects and circumstellar 



emission. 



Wang fc Henning I (120061 ) suggested that o ne of these objects, HH 914, ~ 24" to 



the east of CED 112 IRS 4, is driven by it (see Fig.3 of IWang fc Henning ||200d , where CED 
112 IRS 4 is indicated as Sz 32). Recent ATCA data at 1 6mm show a large con tribution to 
the 16mm flux of CED 112 IRS 4 from free-free emission (jLommen et al. II2009I ). Our HST 
images (Fig. [8]) clearly show that the source is much brighter in the three narrow-band line 
filters than in the broad V-band. In particular, the narrow band images consistently show a 



featur e protruding eastward, which we interpret as the HH 914 object of I Wang fc Henning 



(120061 ) . The Ha image is especially remarkable, as in this filter HH 914 is most clearly 
detached, showing as a second peak ~0 / .'30 (47 AU) to the East about 7.6 times fainter than 
the primary one. The elongation in the V-band to the northeast is most probably due to 
CTE losses. The SED of this source is shown in Section [5j 
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4-0.7. WW Cha (#21) 



WW Cha, a K5 T Tauri star in the vicinity of the previous source CE D 112 IRS4 



( Schegerer et al. 1120061 ). is thought to drive the highly collimated jets HH 915. iBally et al 



(120061 ) report the presence of a giant bow shock, HH 931, about 13' further aw ay from WW 
Cha at P.A.~135°, in the same direction of HH 915. Wang fc Henning] (2006) sugg est that 
two near-infrared if 2 emission knots, A and D, detected by iGomez et al. ( 2004 ) on the 
opposite side of WW Cha may represent the counter-jet of HH 915. There is finally a faint 
chain of Ha knots which links the brightest part of the bow shock to the southern side of the 
reflection nebula illuminated by WW Cha. Studies of the source variability at millimeter 
wavelengths indicate that the 16 mm flux is dominated by cm-size pebbles' emission, which 
makes WW Cha the second star known t o have a protoplanetary disk containing grains 
of suche a large size (jLommen et al. II2009I ). Our images (Fig. [9]) do not show evidence of 



the rich HH phenomenology associated to this source. Hovewer, they show the bright core 
extended in the [SII] filter, with a FWHM contour elongated in the SE-NW direction. If 
the [SII] elongation traces shock emission, then the HH 915 objects would be co-aligned and 
part of the same jet system. The SED of this source is shown in Section [51 



4.0.8. ESO-Ha 569 (#29) 



ESO-Ha 569, spect ral type M2.5, is a high ly variable faint objec t with signatures of bot h 



accretion and outflow (jComeron et al. II2004I ). Also accordin g to IComeron et al. I ( 120041 ) 



the non-detection of this source in the mid-IR by ISOCAM (IPersi et al. 1 120001 ) indicates 
that the amount of warm dust associated with this object is very small, meaning that its 
faintness cannot be attributed to occultation by circu mstellar materi al. On the other hand, 
th e more recent detection of th is source with Spitzer (ILuhman 1 120051 ) and the X-ray image s 
of iFeigelson fc Lawson I ( 120041 ) support the presence of an edge-on disk ( ILuhman 1 120071 ). 
This source is located ~22" to the northeast o f HH 919, an Herb ig-Haro object that could 
be driven by a jet originated by ESO-Ha 569 (IBally et al. 1 120061 ) . In our V-band and Ha 
images, ESO-Ha 569 (Fig. [TO]) is just above our detection limit, but it appears extended in 
the V-band in the NW-SE direction. We do not detect HH 919, nor its associated bipolar 
jet. 
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4-0.9. ESO-Ha 574 (#33) 



ESO-Hct 574 was discovered by IComeron et al. I (120041 ) as a very faint source with the 
characteristic colors of a lightly reddened late-type star. It has a rich emis sion-line spectrum 
domin ated by the forbidden lines typically associated with stellar outflows. IComeron fc Reipurth 
( 120061 ) observed ESO-Ha 574 in the [SII] lines and with the low-resolution spectrograph at 
the ESO-VLT. They detected a well-developed jet (HH 872) protruding from the source and 
discussed the physical proprieties of the emitting source on the basis of the spectrum sam- 
pled at the base of the jet. Recently, iBacciotti et al. I (120111 ) have attributed the unusual low 
luminosity of the source in the continuum to the presence of an edge-on disk. Accounting 
for the disk obscuration on the luminosity of the accretion tracers, they estimate a mass 
accretion rate of ~ 1.7 x 10~ 8 M yr _1 . They also independently derive a mass outflow rate 
in the jet knots of about 1.5 x 10~ 9 , leading to a mass ejection/accretion ratio over the two 
lobes of ~ 0.3. This is in the range expected for magneto-centrifugal jet launch (ICabrit 



20091 ). 



Our WFPC2 images of ESO-Ho; 574 (Fig. [TT]) resolve this source into a nearly edge- 
on disk. In the V-band the disk is bright and extends at P.A.~135°, reaching a length of 
~0'/6 and a thickness, measured at the center of the disk, of ~0'/4. These values correspond 
to 96 AU and 64 AU, respectively, assuming a distance of 160 parsec. The disk and a 
faint trace of a jet are also detectable in the [OI] filter, whereas in the [SII] filter the jet is 
clearly visible, perpendicular to the disk and extended toward t he no rtheast direction. We 
resolved knot A, previously detected by IComeron &: Reipurth I ( 120061 ). in three knots: the 
brightest one, knot A3, stretches ~0'/9 (144 AU) from the center of the disk. The second 
knot, knot A2, is (176 AU) a way from the disk. The third one, knot Al (previously 



identified by IBacciotti et al. 1 1201 lh . is visible at a distance of ~2'/3 (368 AU) from the disk 



and rather than being well collimated as knot A3 it appears bow-s haped. Our [SII] image 



does not show evidence of the southewestern counter-jet resolved by IBacciotti et al. I (120 111 ) 



in their position-velocity diagram, probably because of the relatively low sensitivity of the 
WFPC2/PC to diffuse structures. The SED of this source is shown in Section 



5. Model fitting of the Spectral Energy Distributions 



In order to estimate the main parameters of the star+disk systems combining the IR 
observations (mostly from Spitzer ) with the constraints p osed by our HST imaging survey 
we have used the online tool from iRobitaille et al. I (|2006[ ) to fit the Spectral Energy Distri- 
butions (SEDs) of a number of sources. The tool, based on a grid of 20,000 synthetic models, 
aims at reproducing the SEDs of disks around stars with masses between 0.1 and 5OM and 
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ages between 10 3 and 10 7 years. 



Consistently with these limitations, we have fitted only the SEDs of sources having 
spectral class earlier than M6. We also limited ourselves to sources detected by Spitzer in at 
least one IRAC band, with some flux information at 24/zm or longer wavelengths, typically in 
the sub-mm, as these data points provide a critical constrain to the IR excess from the disk. 
We used the photometric data reported in Tables 4 and 5, plus our own F547M photometry 
presented in Table 3. We excluded from the fit the sources #4, #7 and #25 since these 
are confirmed tight binaries (see Section 4.1 and lLafreniere et al. I (120081 )). Sources #17, 
#27 and #29 have a cluste r member within the beam size of the sub-mm observations 
( IKraus fc Hillenbrand 1 120071 ) and therefore those data points were neglected. 



The SED model fitting tool from iRobitaille et al. I (120061 ) aims at reproducing the pho- 



tometric fluxes by varying a set of 16 parameters, of which the extinction (Ay) and the 
distance (d) vary in ranges defined by the user. The remaining 14 "free" parameters charac- 
terize each model, whose SED is computed at 10 different viewing angles. For each angle of 
each SED, the fitting tool calculates a x 2 to estimate the goodness of the fit. In practice, the 
free parameters are typically constrained by the fluxes at certain wavelengths. For example, 
the mass or the temperature of the star are characterized by the optical region of the SED, 
the inner radius or the flaring angle of the disk by the mid-far IR excess, while to constrain 
the disk mass one needs fluxes at A>100/mi. 

Given the lack of far-IR data, and the uncertainties due to the non-simultaneity of 
the observations, our problem is under-constrained and the SED model fitting tool returns 
several possible solutions. To clean-up the ensemble, for each source we considered only the 
solutions corresponding to values of T e ff and Ay falling within ±100K and ±lmag from 
the values given in Table 4. This is appropriate, considering that the spectral types and 
absolute temperature scale of low-mass stars have typical uncertainty of one sub-class. We 
also constrain the distance to be in the range between 130pc and 190pc, an artifice to allow 



for some variability of the source luminosity (see e.g. iMorales-Calderon et al. 1120111 ). As 



the average of the calculated distances remains close to 160 pc, this does not introduce 
systematic trends in the derived luminosities. We finally require that the difference between 
the x 2 an d the best x 2 of each sou rce must be smaller than 3. This latter criterion, also 
adopted by IRobitaille et al. I ( 120071 ) to fit low-mass young stellar object (YSO) SEDs, has 
been necessary because the model grid is too sparse to reliably search for the true minima 
of the x 2 hypersurface. 

Since the minima of the x 2 hypersurface are difficult to resolve, we could not reliably 
identify the best fit model that formally represe nts the physics of the sou rce, with its asso- 
ciated confidence interval. However, as shown by I Grave fc Kumar I ( 120091 ). the distributions 
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of fundamental parameters such as the stellar mass, age and total luminosity tend to show 
peaks independently on the finesse of the model grid. Following these authors, we have used 
these distribution (constrained by our knowledge of the stellar parameters) to derive the 
"best fit" parameters with their uncertainties, computing a weighted mean and a weighted 
standard deviation for each parameter and each source. For the weights we used the inverse 
of the x 2 returned by the SED model fitting tool. The mean and the standard deviation were 
computed on a logarithmic scale because the parameters in the grid are usually uniformly 
sampled in logarithmic scale. The resulting distributions of the stellar age and mass show 
that when a large number of models can be found, they tend to produce a rather narrow peak 
in the distribution of stellar parameters. Only for three sources, #18, #20 and #29, the 
number of models passing our selection criteria is small and the corresponding histograms 
do not show a peak. Following this approach we were able to fit 19 SEDs, shown in Fig. [T2J 
The physical parameters derived from the SED fitting, together with their uncertainties, are 
listed in Tables [6] and CD 



5.1. Parameters derived from SED fitting 



Table [6] provides the parameters related to the central source: the entry number we 
assigned to each star in Table 2 (column 1); the \ 2 per data-point of the best fit (column 
2); the number of fits that satisfy our selection criteria on T e ff, Ay and x 2 (column 3); 
and finally, in the last four columns, the bolometric luminosity, age, mass and radius of the 
central source. Most of the sources have sub-solar mass down to 0.12 M , a value close 
to the low-mass limit of the fitting tool. The median age, 2.5 Myr, is in agreement with 
previous estimates, while the spread between is 0.5 Myr and 5 Myr, with a few outliers. 
In gener al, this spread is in line w ith what typically found using isochronal timing of PMS 
clusters ([Hillenbrand et al. II2008I ). It is known that seve ral factors may af f ect the estimate 



of th e absolute stellar luminosity ofPMS stars (see e.g. lHartmann I l200ll ; iReggiani et al. 



201 ll ). and in particular edge-on disks can explain the most extreme cases of sub-luminous 
sources. If the age spread is real, star for mation in Cham aele on I proceeded r a ther s lowly 
and may still be ongoing, as suggested by iLuhman I ( 120071 ) and lBelloche et al. I (120 111) . 



Table [7] refers to the main disk parameters. For each source we list the e ntry number (col- 
umn 1); the infrared spectral index a 2 _ 2 4, defined as a = dlog(\F\) / dlog A ( jLada fc Wilking 
19841 ). calculated between 2.2/z m and 24/xm and dereddened using the extinction from Table 



4 and the reddening law from iFlaherty et al. I ( 120071 ) (column 2); the mass of the disks 



either derived by the fitting tool (column 3) or directly estimat ed from the data at 87 0/xm 
(column 4); the dust sublimation radius R su b = -R*(T*/T sn fe) 21 (IRobitaille et al. 1120061 ) for 
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a dust sublimation temperature, T su b = 1600 K (column 65; the inclination of the disk to 
the line of sight (column 6). The SED fitting tool also estimates disk scale height factors 
z factor ~ 0.8 and disk flaring parameters /3 ~ 1.0. These values turn out to be similar for all 
sources within about 10%. 

In Table [7] we provide the model disk mass Mf ED only for sources with some flux 
measure at A>100/im, which allows constraining the amount of dust in the outer disk regions 
where most of the mass resides. These estimates can be compared with the the disk mass 
estimated from the flux density of sub- mm continuum dust emission (column 4), calculated 
assuming optically thin emission at long-wavelengths. The flux density in the sub-mm, 
-fsub-mm, can be converted into an estimate of the disk mass through the relation Mdi s k — 
d 2 F suh ^ mm / (k u B u (T c )), where d is the distance (160 pc), k v is the total (gas + dust) mass 
opacity, B U (T C ) is the Planck function at the characterist ic temperature T c of th e emitting 
dust. For the mass opacity we adopted the opacity law of iBeckwith et al. I ( 1l990l ). i.e. k v = 
0.1 x (z//10 12 Hz) cm 2 /g, whereas we considered a val ue of 20 K for the characteristic outer 
disk temperature (see e.g. lAndrews fc Williams 1120051 ). Unfortunately, the sub-mm data can 
be reliably used to derive the disk mass only for 3 sources, as in th e other cases we have 



either a no n-detection or a cluster member, according to the lists of iKraus fc Hillenbrand 



(120071 ) and lLafreniere et al. I (120081 ). close enough to potentially affect the measured sub-mm 
flux. The values obtained in those cases can be regarded as upper limits to the real disk 
mass, and are therefore reported in square brackets. 



Table 6. Stellar parameters obtained from the fit 



N. x 2 best fit Accepted Fits 


Lbol 


Age 


M* 


R* 






[log yr] 


[M Q ] 


[Re] 



1 


1.1 


242 


0.09 


± 


0.03 


6.81 


± 


0.13 


28 ± 03 


0.81 


+_ 


0.12 


5 


0.087 


818 


0.10 


± 


0.05 


6.57 


± 


0.13 


19 ± 03 


1.0 


± 


0.2 


9 


0.62 


768 


0.18 


± 


0.06 


6.49 


+_ 


0.14 


31+ 03 


1.2 


± 


0.2 


10 


4.3 


12 


0.2 


± 


0.3 


5.9 


± 


0.6 


17 ± 02 


1.5 


± 


0.8 


11 


0.055 


1524 


0.8 


± 


0.3 


6.5 


± 


0.2 


93 ± 07 


1.7 


± 


0.3 


16 


0.46 


1584 


0.08 


± 


0.03 


6.4 


± 


0.3 


12 ± 02 


1.0 


± 


0.2 


17 


0.80 


226 


0.28 


± 


0.13 


6.4 


± 


0.2 


30 ± 03 

V 7 • 'J V 7 — 1 — V 7 • V 7 'J 


1.4 


± 


0.4 


18 


0.33 


6 


0.6 


± 


0.4 


6.0 


± 


0.3 


31 ± 02 


2.1 


± 


0.7 


19 


3.4 


33 


25 


± 


14 


5.2 


± 


0.2 


1.7 ± 0.4 


9 


± 


2 


20 


0.74 


6 


1.2 


± 


0.3 


5.8 


± 


0.2 


0.40 ± 0.03 


2.8 


± 


0.5 


21 


0.44 


90 


20 


± 


12 


5.3 


± 


0.3 


1.6 ± 0.4 


7 


± 


3 


22 


0.27 


54 


1.2 


± 


0.4 


6.4 


± 


0.3 


0.69 ± 0.08 


1.6 


± 


0.4 


23 


0.042 


1070 


0.39 


± 


0.16 


6.4 


± 


0.2 


0.47 ± 0.03 


1.4 


± 


0.3 


24 


0.41 


1640 


0.08 


± 


0.04 


6.5 


± 


0.2 


0.13 ± 0.02 


1.0 


± 


0.2 


26 


0.33 


679 


0.06 


± 


0.06 


6.5 


± 


0.4 


0.12 ± 0.02 


0.8 


± 


0.4 


27 


0.73 


126 


0.5 


± 


0.3 


6.2 


± 


0.2 


0.39 ± 0.03 


1.8 


± 


0.6 


29 


2.0 


5 


0.8 


± 


0.4 


5.8 


± 


0.2 


0.33 ± 0.03 


2.5 


± 


0.6 


30 


0.076 


1602 


0.10 


± 


0.05 


6.59 


± 


0.15 


0.23 ± 0.03 


1.0 


± 


0.2 


32 


0.030 


614 


0.22 


± 


0.09 


6.38 


± 


0.13 


0.32 ± 0.02 


1.3 


± 


0.2 
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NACO - KS - Feb 2006 



C1 




C2 



NACO — H — Feb 2008 



PC - F675N - Apr 2009 



CI 




CENTER: RA 1 1 04 06.65 OEC -76 27 21.5 



Fig. 5. — CT Cha and the probable companions "CI" and "C2" (NACO's images from the 
ESO's archive and PC image). Our detection upper limits for source "CI" are: F547M: 
m > 22.7; F631N: F x < 7.0 x KT 19 erg cm" 2 s" 1 A" 1 ; F656N: F x < 2.8 x 10" 18 erg cm" 2 
s" 1 A -1 ; F673N: F x < 4.6 x KT 19 erg cm" 2 s" 1 A -1 . 
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Fig. 6.— HH 48 A and B observed in the PC chip. Image centered at R.A .=11:04:23.04, 
Decl.=-77: 18:08.3 (J2000.0). 
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Fig. 7. — Ass Cha T 2-16 observed in the PC chip. Image centered at R.A. =11:04:57. 01, 
Decl. =-77:15:57.3 (J2000.0). The insets at the bottom-left corner show the contours for the 
typical PSF. 
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Fig. 8.— CED 112 IRS 4 observed in the PC chip. Image centered at R.A. = 11:09:53.35, 
Decl. =-76:34:25. 5 (J2000.0). The insets at the bottom-left corner show the contours for the 
typical PSF. 
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Fig. 9. — WW Cha observed in the WF chip. Image centered at R.A. =11:10:00. 03, Decl.=- 
76:34:57.3 (J2000.0). The insets at the bottom-left corner show the contours for the typical 
PSF. 
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ESO-Ha 569 - F656N ESO-Ha 569 - F673N 




Arc Seconds Arc Seconds 



Fig. 10.— ESO-Ha 569 observed in the PC chip. Image centered at R.A. = 11:11:10.70, 
Decl.=-76:41:57.6 (J2000.0). 
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Fig. 11.— ESO-Ha 574 observed in the PC chip. Image centered at R.A. = 11:16:03.08, 
Decl.=-76:24:52.6. 
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Fig. 12. — SEDs of our sample of sources. Data obtained in our HST line-filters have not 
been used for the fit as they may be contaminated by accretion or mass loss. They are 
presented in the plot with black error bars. The data actually used for the fit have red error 
bars. 
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5.2. Analysis of the Best-fit Model Parameters 



In this final section we discuss a number of diagrams useful to address the presence of 
evolutionary trends between the main parameters of our star+disk systems. 

First, as a sanity check, we compare the values of stellar age and mass returned by 
the fitting tool with those obtained from the di rect interpolation of T e *t and L^i in the 



HR diagram using the Siess tracks (ISiess 1 120001 ) , which are the same tracks adopted by 
the fitting tool. Fig. [13] and [14] show that there is strong agreement. In particular, the 
strong correlation between the mass values is expected, since the model fitting assumed 
temperatures within ±100 K from the stellar temperatures and in this temperature range 
the mass is strongly correlated to the temperature. There are, however, four sources (#1, 
#10, #19, #21) for which the fitting tool provides ages and masses which are discrepant 
from those found in the literature. The younger ages we derive are a result of the higher 
bolometric luminosities returned by the fitter. For disk seen at high inclinations, taking into 
account the IR part of the SED allows to recover a non- negligible fraction of the stellar flux 
which would otherwise rem ain unaccounted for b y a applyi n g a si mple reddening correction 
to the optical photometry. iLuhman I (120041 ) and iLuhman I (120071 ). who derived Lboi from I 
and J band magnitudes for almost the entire known population of Chamaeleon I, already 
noticed that sources with unusually low L&„j values, which apparently lie below the Main 
Sequence in the HR diagram, may be highly affected by an underestimate of the dust column 
density. Unfortunately these systems are also the most challenging to model, due to the lack 
of direct information e.g. on the disk flaring angle and on the dust properties at the disk 
surface. With the exception of source #10, whose age (i.e. luminosity) derived from the fit 
is highly uncertain, the other three sources are distributed at the two extremes of the age 
range: source #1 is the most luminous while source #19 and #21 are the faintest ones. Our 
HST images show that these last two objects are associated with diffuse emission, consistent 
with sources seen mainly in scattered light. 

Table [7] shows that 14 out of 19 disks have tilt angle within 30 degrees from edge-on, 
about twice the number one would expect if disks are oriented randomly. While this may 
be partially due to the sparse sampling of this parameter provided by the theoretical model, 
we must remark that the model may actually be biased toward larger angles, as it may 
not sample disks that are enough settled. The model adopts a fully flared disk (in vertical 
hydrostatic equilibrium) and mimics the dust settling by multiplying the sca le height at 
the dust sublimation radius by a "zfactor" . Figure 6 in iRobitaille et al. I ( 120061 ) shows that 
the zfactor is a function of the disk outer r a dius, becoming equal or less than 0.5 only for 
relatively large disks. However, ISzucs et al I (120101 ) notice that the average SED of T-Tauri 
stars in the Chamaeleon-I region already requires a reduction of a factor of 2 in the disk scale 
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Table 7. Disk parameters obtained from the fit 



N. 


«2-24 


m sed 


j^j-sub—mma 




Inclination 








[Mq] 


[AU] 


[cleg] 


f 


-0.52 ± 0.12 


— 





0.02 


81.4 


5 


-1.14 ± 0.12 


— 


<0.13 


0.02 


75.5 


9 


-1.35 ± 0.12 


— 


— 


0.03 


75.5 


fO 


-0.60 ± 0.12 




0.003 


0.03 


87.1 


ff 


-1.27 ± 0.27 


— 


— 


0.06 


63.3 


16 


-2.45 ± 0.25 


— 


— 


0.016 


87.1 


17 


-0.53 ± 0.14 




[0.12] 


0.03 


18.2 


18 


-0.29 ± 0.17 






0.05 


87.1 


19 




0.001 ± 0.004 


0.005 


0.3 


75.5 


20 


0.20 ± 0.17 




< 0.03 


0.8 


87.1 


21 




0.02 ± 0.08 


0.12 


0.3 


31.8 


22 


-0.66 ± 0.12 




[< 0.03] 


0.05 


41.4 


23 


-1.14 ± 0.12 






0.04 


69.5 


24 


-2.58 ± 0.17 






0.016 


81.4 


26 


-0.95 ± 0.12 






0.016 


75.5 


27 


-0.59 ± 0.12 




[0.003] 


0.05 


18.2 


29 


-0.08 ± 0.19 




[0.005] 


0.05 


87.1 


30 


-1.32 ± 0.14 






0.02 


75.5 


32 


-2.69 ± 0.13 






0.03 


56.6 



a In this column we report the disk mass derived from the data at 
870 ^m (see Section 2.3). The values relative to sources with a visual 
companion within the beam are given in brackets and can be considered 
as upper limits to the real disk masses. 
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Fig. 13. — Comparison between the a ges estimated with Luhm an luminosities + Siess tracks 
and the SED model fitting tool from iRobitaille et al. I ( 120061 ). 
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Fig. 14. — Comparison between the stellar mass es estimated with Luhm an luminosities + 
Siess tracks and the SED model fitting tool from IRobitaille et al. I ( 120061 ). 



height; similar r esults are also repo rted by the Spitzer Infrared Spectrograph GTO team for 
Taurus sources (IFurlan et al.ll201ll ). This suggests that the grid of Robitaille et al. may not 
properly sample the typical settling of TTau disks. In conclusion, the returned values for 
the tilt angle are only indicative, and therefore we use them to flag outliers without drawing 
any further conclusion. 



To further investigate the possible correlations between the disk structure and the de- 
rived stellar properties, we compare the infrared spectral index 02-24, a tracer of the amount 
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-3 



-2 



-1 

^2.2-24 



Fig. 1 5. — Difference between the values of L^i derived in this work and those from iLuhman 
( 120071 ) versus the infrared spectral index 02-24 that is a tracer of the structure of the inner 
disk, being mostly dependent on the flaring angle and on the presence of an inner rim. 



of dust in the warmest layers of the circumstellar disk, to the main parameters of the central 
sources, like luminosity and mass. Fig. [T5l s hows the difference between the values of L^oi 



derived by the fitting tool and by ILuhman I ( 120071 ) against the spectral index 02-24- Three 



sources with discrepant luminosity stand out in the plot: #18, #20 and #29. T hese three 



sources, with nearly flat IR SEDs, were positioned according to ILuhman I (120071 ) below the 
Zero- Age-Main-Sequence and therefore, having no estimate for their mass and age, did not 
appear in Figures [13] and HU The fitting tool returns for them a small number of successes, 
confirming that these systems have extreme properties. In particular, these are the sources 
with the highest tilt angle, seen nearly edge-on. Our derived ages and mass are now in-line 
with the main population (see Figures [161 and [T71 below). Sources #19 and #21 do not 
appear in this plot as they have no detection at 24 /im. 

Fig. [Tj2 also shows that for the other sources the difference between the estimated lu- 
minosities appears to increase with the spectral index: there is a correlation between the 
amount of warm dust seen in the outer disk layers and the apparently low bolometric lumi- 
nosity derived from dereddening the I and J-band magnitudes. In general, this plot seems 
to confirm an underestimate of the dust column density toward the edge-on sources. 

Fig. [T6l shows that the index 02-24 decreases with time for nearly all sources (#19, and 
#21, younger than 0.2 Myr and without a 24 /im measure, are not plotted), indicating that 
the disks become flatter with age. This is most probably related to the gradual settling of 
the dust grains on the disk plane. Our linear fit to the main distribution is provided by 
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the formula shown in the figure. iLuhman et al. I ( 120081 ) set 02-24 = —2.2 as the lower limit 
for disk dissipation. According to our expression, this value is reached at a time ~ 10 7 yr, 
compatible with the scenario that disks around mo st young solar analog stars cle an out their 
small dust grains within 1 AU in w 10 Myr or less (IPascucci fc Tachibana II2010I ). The three 
sources that strongly deviate from the main trend, #16, #24 and #32, are those exhibiting 
negligible IR excess. These outliers seem to have dissipated their inner disk relatively quickly, 
less than ~ 3 Myr. Various effects may lead to rapid disk dissipation, all intimately linked to 
the physical or chemical properties of the environment, like e.g. higher UV or X ray flux from 
the central star, ti dal forces due to the presence of a close companion or of a giant planet 
(ICieza et al. 1120091). or low disk metallicity which can increment the dust photoevaporation 
rate (jErcolano et al. 1120101 ). A larger statistical sample is needed to quantify the frequency 
of premature versus delayed disk dissipation. 

A plot of the «2-24 index against stellar mass (Fig. [I?]) , does not show any clear corre- 
lation. However, two of the three sources at the bottom of the plot, without an inner disk, 
have stellar mass close to our lower limit. For this sources inner disk dissipation may be 
more likely driven by tidal forces or rapid formation of a giant planet rather than photoevap- 
oration, which is expected to be dominant in more massive sources. ILuhman et al. I (120081 ) 
noticed that Chamaeleon I, unlike other star forming regions, contains a significant fraction 
of low-mass stars with inner-disk lifetimes shorter than those of more massive stars. 
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Fig. 16. — Infrared spectral index 0:2-24 as a function of the age. Black dots represent the 
sources used to fit the main trend (dashed line) and red dots represent the diskless sources. 

In what concerns the mass accretion rates derived from our HST Ha photometry, sev- 
eral observations have shown a correlation between the mass accretion rate and the age, 
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Fig. 17. — Infrared spectral index «2-24 in function of the mass. Red dots represent the 
diskless sources. 



mass or the 1 


-1 spectral index of YSOs ( 


Muzerolle et al. 


2003, 


2005; 


Mohantv et al. 


2005; 


Natta et al. 1 


2006 


: I 


Alexander & Armitaee 1 


2006 




Gatti et al. 


2008 


: ISicilia-Aeuilar et al. 


2010 




Manara et al. 


2012). Based on the values reported in Table IB1 and the magnitudes 



listed in Table H] we investigate the correlation between these parameters and the mass accre- 
tion rates estimated from our HST images. For sources #18, #20 and #29 we use our new lu- 
minosity estimates (Table [6]) to rederive the mass accretion rates, obtaining logM acc = —8.2, 
-7.9 and -7.38 M /yr _1 , respectively, and use these values instead of those presented in Ta- 
ble El We also discard sources # 7 and #25, as they are confirmed tight binaries (see Section 



4.1 and lLafreniere et al. 1 120081 ) 



In Fig. [18] we plot the mass accretion rate as a function of stellar age. The largest 
accretion rates are found for the two youngest stars, while the majority of older sources 
show a spread of about 2 orders of mangnitude. These characteristics, a general decrease 
of the mass accretion rate vs. time associated with a large scatter at any given age, cannot 
be explained by a any reasonable systematic overestimate of the stellar luminosity (which 
would make the stars younger while enhancing the estimated accretion lu minosity) and are 
consi stent with what typically found in other star forming regions (see e.g. lHartmann et al. 
19981 ). 



In Fig. [19] we show the mass accretion rate as a function of the stellar mass. The plot 
shows a gradual rise of the mass accretion rate with the mass of the centr al source. The 
dashed line sh ows the M aca oc M, 2 s caling relation repo rted by various authors ( jMuzerolle et al. 



20031 . l2005t iMohantv et al. 1 120051 ; iNatta et al. 1 120061 ) and discussed e.g. in the context of 
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Fig. 18. — Mass accretion rate plotted as a function of the stellar age. The age values are 
those extrapolated from the SED fitting tool and listed in Table O 



accelerated disk cle aring by Oarke fc Pringle I (20061 ) or of systematic differences in disk 
initial conditions by [Alexander fc Armitage I ( 120061 ). 



age < 1 Myr 

age > 1 Myr ^ 
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Fig. 19. — Mass accretion rate plotted as a function of the stellar mass. The mass of the 
central sources are those extrapolated from the SED fitting tool and listed in Table |6j The 
dots represent the sources with age>l Myr, whereas the younger sources are indicated by a 
triangle. The dashed line represent the relation M acc (xM^. 



Finally, in Fig. [20] we show the relation between the mass accretion rate and the IR 
spectral index 0:2-8, more sensitive to the warmer dust in the inner disk than the 0:2-24 
index. The 02-8 index is a tracer of both disk evolutionary status and tilt angle. While the 



scatter of points in Fig. [20] may be attributed to the tilt angle, the systematic trend suggests 
that more evol ved inner disks tend to have smaller accretion rates. A similar correlation has 
been found by ISicilia-Aguilar et al. I (120 lOh . 




Fig. 20. — Mass accretion rate plotted as a function of the IR spectral index a 2 _s- 



6. Conclusions 

We have used the WPC2 instrument onboard the HST to observe 20 fields centered on 
T Tau stars in the star forming region Chamaeleon I. Our images, obtained in narrow-filters 
centered on the [01] , Ha and [SII] lines, plus a Johnson- V band equivalent filter, allowed us 
to detect 31 previously known T Tauri stars. In this paper we have presented the images 
relative to 10 sources that appear either extended, binaries, or surrounded by a circumstellar 
disk and/or mass outflow. We have complemented our photometry with a compilation of 
optical, IR and sub-mm observations, adding new sub-mm data for three sources, together 
with published values of T e ff, Aj and L hoi . Using our Ha photometry we have estimated 
the mass accretion rates for 28 sources. Using all available data, we have reconstructed the 
optical-IR SED s of 19 sources and der ived a number of disk parameters using the SED-model 



fitting tool of iRobitaille et al. I ( 120071 ). Our main results are: 



1. We resolved 10 Chamaeleon I sources into binaries or diffuse objects with evidence of 
circumstellar material, either disks or jets. 



2. The SED fitting shows that most of the sources have mass between 1.7 M and 0.12 M Q 
and isochronal ages typically ranging between 0.5 Myr and 5 Myr. 
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3. The Lhni d erived from the fitting tool is generally higher than the values reported by 



Luhman I (120071 ). leading to apparently younger ages. The discrepanc y increases with 



the I R spe ctral index Qj^_ ?4 (see Fig. [TBI) . A few sources identified by iLuhman I (12004 



2007 ) and Luhman &: Muench (j2008l ) as anomalously faint in the near-IR appear to 
be affected by high dust column density and most probably are associated with disks 
seen nearly edge-on. 

4. For 13 sources the IR spectral index «2-24 appears to decrease with time. The rate of 
decay would imply that disk dissipation in Chamaeleon I requires ~ 10 7 yr, a relatively 
long disk lifetime. There are 3 sources that deviate from the main trend, exhibiting an 
unusual absence of IR excess for their apparently young age. 

5. The mass accretion rate decreases with the stellar age, showing a spread of about two 
orders of magnitude at ~ 3 Myr, consistently with what typically found in other star 
forming regions. The mass accretion increases with the stellar mass roughly following 
the same scaling relation, M acc oc Ml found in other PMS clusters. 

6. The fact that both the IR spectral index «2-24 and the mass accretion rates decrease 
with our estimated isochronal time suggests that the age spr ead observed in our s ample 



is real. This is in contrast with the recent suggestion by iJeffries et al. I (120111 ) that 
individual stellar ages from the Hertzprung-Russel diagram are unreliable since, at 
least in the Orion Nebula Cluster, they do not correlate with the presence of disks 
inferred from near-IR excess. There are clearly several factors that may contribute to 
the observed luminosity dispersion, and we have shown that the case of highly tilted 
disks is one of those. A comprehensive analysis of the rich phenomenology associated 
with the presence of accreting circumstellar disks may allow to reveal the intrinsic age 
spread within a cluster. 
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A. SEDs of individual sources 



2MASS J10533978-7712338 (#1): for this source the SED fitting tool is unable to re- 
produce the relatively flat slope of the Spitzer/IRAC data (3.5-8.0/xm). The absence of 
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near-IR excess drives the solution toward a highly tilted system, i.e. a disk with 81.4° 
inclination with r espect to the plane of the sky. The high tilt angle agrees with pre- 
vious suggestions (ILuhman &: Muench 1120081 ) that this source (unresolved in our HST 
images) is probably mainly seen in scattered light. The low bolometric luminosity 
makes this source the oldest one of our sample. 

ISO-Chal 10 (#5): for this source we obtain a good fit, with st e llar m ass is in agreement 



with the value of ~0.18 M report ed by lLafreniere et al. I (120081) . This source has 



20041 ) but remains unresolved 



been indicated as a possible binary (ILopez-Marti et al. 
in our HST observations. 

Ass Cha T 2-16 (#9): also for this source we obtain a a good fit, The stellar mas s 0.19 M ffl 
is smaller than the previous estimate of 0.26 M by lLafreniere et al. I (120081 ). In Sec- 
tion ?? we showed some evidence for a spatially resolved PSF, especially in the [OI] 
line filter. 

ESO Ha-559 (#10): for this source we have a sub-mm detection. The best fit indicates 
a disk seen nearly edge-o n (at 87.1° tilt) with a disk mass ~2-10~ 3 M Q and an age 
~7.9 Myr, consistent with lComeron et al. I (120041 ). 



V* UX Cha (#11): the SED of this source is typical of a transition disk that has almost 
entirely dissipated the inner region. The fit returns a stellar mass of ~0.9 M , the 
same value assigned by Kirk fc Myers I ( 120111 ) . 



ISO-Chal 201 (#16): the SED of this source is compatible with a pure photosphere up 
to 24 /im. The inner disk has been cleared rather rapidly, as the source seems only 
~2.5 My old. This source has been classified as a candidate brown dwarf with spectral 
type M5.75; our fit assigns a mass of ~0.12 M , close to the lower limit of the grid 
values. 

Hn 10E (#17): the best fit for this source provides a marginally acceptable match to the 
flat IR SEP. The derived stellar parameters are consistent with those reported by 



Feigelson fc Lawson I (120041 ) . 



CHSM 15991 (#18): even if the fitting tool provides a small number of acceptable solu- 
tions, the best fit for this source shows good agre ement with the data. W e derive an 



extreme disk inclination, 87.1°, in areement with iLuhman et al. I (120081 ) . Cases like 



this of extre me disk in c linati on make the estimates of the absolute stellar lumiosity 
problematic. ILuhman I (120071 ) estimated for this source Lb Q i ~ 0.0029 L , which puts 
below the Main Sequence it in the H-R diagram. Our best fit returns a much higher 
luminosity, Lfc o z~0.6 L and an age ~1.0 Myr, for a ~0.31 M stellar mass. 
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CED 112 IRS 4 (#19): we obtain a good fit for this source, resolved in our HST images 
and detected at 870/im. The young age, 0.14 Myr, is well compatible with the HST 
images showing a young active source associated with HH 914. The disk inclination, 
75.5°, seems also compatible with the HST morphology. 

ISO-Chal 225 (#20): for this source we find a small number of acceptable models. In 
fact, our best fit poorly reproduces the near-IR photometry. The extreme tilt angle, 
87.1°, would imply that the source, unresol yed in our HST images, is mainly seen in 
scattered light. This may well be the iLuhman I ( 120071 ) estimates L;, o ;~0.013 L© 

putting the star below the mean sequence. We derive a much higher value, L feo ;~1.2 L , 
thus an age of ~0.59 Myr for a mass ~0.40 M . 

V* WW Cha (#21): this source, marginally resolved in our HST images, has two detec- 
tions in the mm region. SED fitting provides a disk mass of ~0.02 M Q and an age 
of ~2-10 5 yr. As for CED 112 IRS 4, high dust extinction may explain the difference 
between a stellar mass ~1.6 M returned from the fit and the estimate of ~0.7 M Q by 



Lafreniere et al. I (120081 ) 



ISO-Chal 232 (#22): we obtain a good fit for this source, which shows strong IR excess 



and mass los s, being associated with objects HH 917 (IBally et al. 1 120061 ) . HH 912 



and HH 916 (jWang fc Henning 1 120061 ) . The stellar mass r eturnd by the fitt e r, 0.6 9 ± 
0.09 M Q , is slightly higher than the 0.55 M Q estimated by lLafreniere et al. I (120081 ). 



Cha T 2-45a (#23) : we obtain a good SED fit. The derived stellar mass, 0.4 7±0.03 M e 



is in agreement with the estimate (~0.51 M ) of lLafreniere et al. I (120081 ). 



0; 



Hn 12W (#24): The SED shows no evidence of IR excess up to 24/xm. The stellar mass 
0.13± 0.02 M is again in agreement with the estimate (~0.15 M ) of lLafreniere et al. 



ISO-Chal 252 (#26): we obtain a genarally good fit except for the IRAC 8 micron data 

point. The stellar mass, 0.12±0.02 Mm is in agreement with the estimates of lLafreniere et al. 
(120081 ) and iMuzerolle et al. I ( 120051 ). lying out the borderline between stars and brown 
dwarfs. 

HBC 584 (#27): for this source with incompled Spitzer/IRAC coverage and strong 02-24 
index, the best fit indicates a stellar mass 0.39 ± 0.03 M , in agreement with the 



estimated by lLafreniere et al. I ( 120081 ) of 0.35M, 



ESO Hct-569 (#29): a problematic fit for this faint source at optical wavelenghts. The 
SED does not show a strong excess at A < 10 /xm, but th e 24/xm data point is remark - 
ably high. For this source, non detected in the X-rays by iFeigelson &: Lawson I (120041 ). 
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an extinction A^>,60, possibly due to an edge-on circumstellar disk has been suggested 
by iLuhman I (120071 ). The steep rise of the flux at A < 10 /mi confirms that the low 
optical luminosity can be attributed to high dust column density Our best fit supports 
this scenario, with a ~87.1° disk inclination. The L^~0.003 L derived by ILuhman 



( 120071 ) places this source below the Main Sequence in the H-R diagram, but integrating 
the emission reprocessed at longer wavelengths our fit provides L& / ~ 0.1 L Q , with an 
age ~ 0.63 Myr and a mass ~ 0.33 M . 

Hn 21 W (#30): we obtain a good fit to the SEP with a stellar ma s s 0.2 3 ± 0.03 M in 
agreement with the 0.20 M & value reported by lLafreniere et al. I (120081 ) . 



BYB 53 (#32): a class HI SEP. The age ret urned by the fit is consistent with the 2 Myr 
estimated by I Gomez fc Mardones I (120031 ). indicating that the source has depleted its 
disk very quickly. 



REFERENCES 

Alexander, R. P., Armitage, P. J. 2006, ApJ, 939L, 83 

Allard, F., Homeier, P., Freytag, B. 2010, 2010arXivl011.5405A 

Allen, L. E., et al. 2004, ApJS, 154, 363 

Andrews, S. M. & Williams, J. P. 2005, ApJ, 631, 1134 

Andrews, Sean M., Wilner, P. J., Espaillat, C, Hughes, A. M., Pullemond, K., McClure, 
M. K., Qi, C. & Brown, J. M. 2011, ApJ, 732, 42 

Apai, P., Pascucci, I., Bouwman, J., Natta, A., Henning, T., Pullemond, C. P. 2005, Science, 
310, 834 

Apai, P., Lauretta, P. S. 2010, Protoplanetary Pust: Astrophysical and Cosmochemical 
Perspectives, Cambridge University Press 

Arce, H., Shepherd, P., Gueth, F., Lee, C. F., Bachiller, R., Rosen & A., Beuther, H. 2007 
in Protostars and Planets V, University of Arizona Press, Tucson, 951, 245 

Bacciotti, F., Whelan, E. T., Alcala, J. M., Nisini, B., Podio, L., Randich, S., Stelzer, B. & 
Cupani, G. 2011, ApJ, 737, 26 

Baggett, S., et al. 2002, in HST WFPC2 Pata Handbook, v. 4.0, ed. B. Mobasher, Baltimore, 
STScI 



-52 - 



Bally, J., O'Dell, C. R., McCaughrean, M. J. 2000, AJ, 119, 2919 

Bally, J, Walawender, J, Luhman, K. L., Fazio, G. 2006 AJ, 132, 1923 

Bally, J., Reipurth, B., & Davis, C. J. 2007, in Protostars and Planets IV, Univ. of Arizona 
Press, 215 

Beckwith, S. V. W., Sargent, A. L, Chini, R. S., Guesten, R. 1990, AJ, 99, 924 

Belloche, A., Schuller, F., Parise, B., Andre, P., Hatchell, J., Jorgensen, J. K., Bontemps, 
S., Weiss, A., Menten, K. M. & Muders, D. 2011, A&A, 527, 145 

Biretta, J. & McMaster, M., ed 2008, Wide Field and Planetary Camera 2 Instrument Hand- 
book, Version 10.0, (Baltimore: Space Telescope Science Institute) 

Blum, J. & Wurm, G. 2008, ARA&A, 46, 21 

Bouvier, J., Alencar, S. H. P., Harries, T. J., Johns-Krull, C. M., Romanova, M. M. 2007, 
in Protostars and Planets IV, Univ. of Arizona Press, 479 

Cabrit, S. 2009, in Protostellar Jets in Context , by K. Tsinganos, T. Ray, M. Stute. Astro- 
physics and Space Science Proceedings Series, Berlin, Springer p. 247 

Calvet, N., D'Alessio, P., Hartmann, L., et al. 2002, ApJ, 568, 1008 

Cambresy, L. et al. 1997, A&A, 324, L5 

Cambresy, L. et al. 1998, A&A, 383, 892 

Carpenter J.M. , Hillenbrand L.A., Skrutskie M.F., Meyer M.R. 2002, AJ, 124, 1001 
Cieza, L. A. et al. 2009, ApJ, 696, L84 

Clarke, C. J., Gendrin, A., Sotomayor, M. 2001, MNRAS, 328, 485 
Clarke, C. J. & Pringle, J. E. 2006, MNRAS, 370, L10 

Colgan, S. W. J., Schultz, A. S. B., Kaufman , M. J., Erickson & E. F., Hollenbach, D. J. 
2007, ApJ, 671, 536 

Comeron, F., Reipurth, B., Henry, A., Fernandez, M. 2004, A&A, 417, 583 

Comeron, F., Reipurth, B. 2006 A&A, 458, L21 

Da Rio, N., Robberto, M., & Soderblom, D. R. 2010, ApJ, 722, 1092 



-53 - 



De Marchi, G., Panagia, N., Romaniello, M. 2010, ApJ, 715 ,1 

Dolphin, A.E., 2009, "A revised Characterization of the WFPC2 CTE Loss", Astronomical 
Society of the Pacific, 121:655-667 

Dominik, C, Blum, J., Cuzzi, J. N. & Wurm, G. 2007 in Protostars and Planets IV, University 
of Arizona Press, Tucson, 951, 783 

Ercolano, B., Clarke, C. J. 2010, MNRAS, 402, 2735 

Espaillat et al. 2008, ApJ, 682, 125 

Feigelson, E. D., Lawson, W. A. 2004, ApJ, 614, 267 

Flaherty, K. M.,Pipher, J. L., Megeath, S. T., Winston, E. M., Gutermuth, R. A., Muzerolle, 
J., Allen, L. E. & Fazio, G. G. 2007, ApJ, 663, 1069 

Fruchter, A. and Sosey, M. et al. 2009, The MultiDrizzle Handbook, Version 3.0, (Baltimore: 
Space Telescope Science Institute) 

Furlan, E., Luhman, K. L., Espaillat, C, D'Alessio, P., Adame, L., Manoj, P., Kim, K. H., 
Watson, Dan M., Forrest, W. J., McClure M. K., Calvet, N., Sargent, B. A., Green, 
J. D., Fischer, W. J., 2011, ApJS, 195, 3 

Gatti, T., Natta, A., Randich, S., Testi, L., Sacco, G. 2008, A&A, 481, 423 

Gauvin, L. S., Strom, K. M. 1992, ApJ, 385, 217 

Gomez, M., Mardes, D. 2003, A J, 125, 2134 

Gomez, M., Persi, P., Marenzi, A. R., Roth, M., Tapia, M. 2004, A&A, 423,629 

Grave, J. M. C. & Kumar, M. S. N. 2009, A&A, 498, 147 

Gullbring, E., Hartmann, L., Briceno, C, & Calvet, N. 1998, ApJ, 492, 323 

Gutermuth, R. A., Megeath, S. T., Muzerolle, J., Allen, L. E., Pipher, J. L., Myers, P. C, 
Fazio, G. G. 2004, ApJS, 154, 374 

Haisch, K. E., Jr., Lada, E. A. & Lada, C. J. 2001, ApJ, 553, L153 

Haisch, K. E., Jr., Greene, T. P., Barsony, M., & Stahler, S. W. 2001, ApJ, 553, L153 

Hartmann, L. 2001, AJ, 121, 1030 



-54 - 

Hartmann, L., Calvet, N., Gullbring, E., D'Alessio, P. 1998, ApJ, 495, 385 

Henize, K. G., Mendoza, E. E. 1973, ApJ, 180, 115 

Henning, T., Pfau, W., Zinnecker, H. & Prusti, T. 1993, A&A, 276, 129 

Hillenbrand, L. A., Bauermeister, A., White, R. J. 2008, 14th Cambridge Workshop on Cool 
Stars, Stellar Systems, and the Sun ASP Conference Series, Vol. 384, 200 

Jeffries, R. D., Littlefair, S. P., Naylor, T., & Mayne, N. J. 2011, MNRAS, 418, 1948 

Kessler-Silacci, J. et al. 2006, ApJ, 639, 275 

Kim, K. H. et al. 2009, ApJ, 700, 1017 

Kirk, H. & Myers, P. C. 2011, ApJ, 727, 64 

Kraus, A. L., Hillenbrand, L. A. 2007, ApJ, 662, 413 

Krumholz, M. R. 2006, ApJ, 641, 45 

Lada, C. J. & Wilking, B. A. 1984, ApJ, 287, 610 

Lafreniere, D. et al. 2008, ApJ, 683, 844 

Lawson, W. A., Feigelson, E. D., and Huenemoerder, D. P. 1996, MNRAS, 280, 1071 

Lommen, D. Maddison, S. T., Wright, C. M., Van Dishoeck, E. F., Wilner, D. J., Bourke, 
T. L., 2009 A&A, 495, 869 

Lopez-Marti, B., Eisloffel, J., Scholz, A., Mundt, R. 2004 A&A, 416, 555 

Luhman, K. L. 2004, ApJ, 602, 816 

Luhman, K. L., 2005, ApJ, 631, L69 

Luhman, K. L. 2007, ApJS, 173,104 

Luhman, K. L. 2008, in Handbook of Star Forming Regions, ASP Conf. Series 
Luhman, K. L., Muench, A. A. 2008, ApJ, 684, 654 

Luhman, K. L., Allen, L. E., Allen, P. R., Gutermuth, R. A., Hartmann, L., Mamajek, E. 
E., Megeath, S. T., Myers, P. C, Fazio, G. G. 2008, ApJ, 675, 1375 

Manara, C. F. et al. 2012, ApJto be submitted 



-55 - 



Masciandri, E., Raga, A 2004, ApJ, 615, 850 
Megeath, S. T., et al. 2004, ApJS, 154, 367 

Meyer, M. R. et al. 2006, Astronomical Society of the Pacific, 118:850, 1690-1710 

Millan-Gabet, R., Malbet, F., Akeson, R., Leinert, C, Monnier, J. & Waters, R. 2007, 
Protostars and Planets V, B. Reipurth, D. Jewitt, and K. Keil (eds.), University of 
Arizona Press, Tucson, 951, 539-554 

Mohanty, S., Jayawardhana, R., Basri, G. 2005, ApJ, 626, 498 

Morales-Calder/'on, M., et al. 2011, ApJ, 733, 50 

Muzerolle , J., Calvet , N., Hartmann , L. 2001, ApJ, 550, 944 

Muzerolle, J., Hillenbrand, L., Calvet, N., BricenS, C, Hartmann, L. 2003, ApJ, 592, 266 
Muzerolle, J. et al. 2004, ApJ, 617, 406 

Muzerolle, J., Luhman, K. L., Briceno, C, Hartmann, L., Calvet, N. 2005, ApJ, 625, 906 

Natta, A., Meyer, M. R., Beckwith, S. V. W. 2000, ApJ, 534:838-845 

Natta, A., Testi, L., Randich, S. 2006, A&A, 452, 245 

O'Dell, C. R. & Wen, Z. 1994, ApJ, 436, 194 

O'Dell, C. R. & Henney, W. J., 2008, AJ, 135, 1566 

Padgett, D. et al. 1999, AJ, 177, 1490 

Pascucci, I., Apai, D., Hardegree-Ullman, E. E., Kim, J. S., Meyer, M. R., Bouwman, J. 
2008, ApJ, 673, 477 

Pascucci, I., & Sterzik, M. 2009, ApJ, 702, 724 

Pascucci, I., & Tachibana S., 2010, in Protoplanetary Dust: Astrophysical and Cosmochem- 
ical Perspectives, eds.: D. Apai, D. S. Lauretta, 2010, Cambridge University Press, 
p. 263 

Persi, P., et al. 2000, A&A, 357, 219 

Podio, L., Bacciotti, F., Nisini, B., Eisoffel, J., Massi, F., Giannini, T., Ray, T. P. 2006, 
A&A, 456, 189 



-56 - 

Reggiani, M., Robberto, M., Da Rio, N., Meyer, M. R., Soderblom, D. R., Ricci, L. 2011, 
A&A, 534, 83 

Reipurth, B., Yu, Ka Chun, Heathcote, S., Bally, J., Rodriguez, L. F. 2000, AJ, 120, 1449 
Richer, J.S. et al. 2000 in Protostars and Planets IV, Univ. of Arizona Press, Tucson, p. 867 
Ricci, L., Robberto, M., Soderblom, D. R. 2008, AJ, 136, 2136 

Ricci, L.,Testi, L., Natta, A., Neri, R., Cabrit, S. & Herczeg, G. J. 2010, A&A, 512, A15 

Robberto, M., Song, J., Mora Carrillo, G., Beckwith, S. V. W., Makidon, R. B., Panagia, 
N. 2004, ApJ, 606, 952 

Robberto, M., Beckwith, S. V. W., Panagia, N., Patel, S. G., Herbst, T. M., Ligori, S., 
Custo, A., Boccacci, P., Bertero, M. 2005, AJ, 129, 1534-1563 

Robitaille, T. P., Whitney, B. A., Indebetouw, R., Wood, K., & Denzmore, P. 2006, ApJS, 
167, 256 

Robitaille, T. P., Whitney, B. A., Indebetouw, R., & Wood, K. 2007, ApJS, 169, 328 
Rydgren, A.E. 1980, AJ, 85, 444 

Schegerer, A., Wolf, S., Voshchinnikov, N. V., Przygodda, F., Kessler-Silacci, J. E. 2006, 
A&A, 456:535-548 

Schmidt, T. O. B., Neuhauser, R., Seifahrt, A., Vogt, N., Bedalov, A., Helling, Ch., Witte, 
S., Hauschildt, P. H. 2008, A&A, 491, 311 

Schwartz, R. D., 1977 ApJ, 35:161-170 

Sicilia-Aguilar, A., Henning, T., & Hartmann, L. W. 2010, ApJ, 710, 597 
Siess, L., Dufour, E. & Forestini, M. 2000, A&A, 358, 593 
Siringo, G., et al. 2009 A&A, 497: 945-962 

Smith, N., Bally, J., Shuping, R. Y., Morris, M., Kassis, M. 2005, AJ, 130, 1763 
Szucs, L., Apai, D., Pascucci, I., Dullemond, C. P. 2010, ApJ, 720, 1668 
Tuthill, P. G., Monnier, J. D. & Danchi, W. C. 2001, Nature, 409, 1012-1014 
Wang, H. & Henning, T. 2006, ApJ, 643: 985 



Whelan, E. T., Ray, T. P., Bacciotti, F., Natta, A., Testi, L., Randich, S. 2005, Nature, 435, 
652 



This preprint was prepared with the A AS IATgX macros v5.2. 



